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Abstract 

TCAT, originally developed for the Space Systems Design Lab at the Georgia 
Institute of Technology is a conceptual design tool capable of integrating aeroheating 
analysis into conceptual reusable launch vehicle design. It provides TPS unit thicknesses 
and acreage percentages based on the geometry of the vehicle and a reference trajectory 
to be used in calculation of the total cost and weight of the vehicle design. 

TCAT has proven to be reasonably accurate at calculating the TPS unit weights 
for in-flight trajectories; however, it does not have the capability of sizing TPS materials 
above cryogenic fuel tanks for ground hold operations. During ground hold operations, 
the vehicle is held for a brief period (generally about 2 hours) during which heat transfer 
from the TPS materials to the cryogenic fuel occurs. If too much heat is extracted from 
the TPS material, the surface temperature may fall below the freezing point of water, 
thereby freezing any condensation that may be present at the surface of the TPS. 
Condensation or ice on the surface of the vehicle is potentially hazardous to the mission 
and can also damage the TPS. It is questionable whether or not the TPS thicknesses 
provided by the aeroheating analysis would be sufficiently thick to insulate the surface of 
the TPS from the heat transfer to the fuel. Therefore, a design tool has been developed 
that is capable of sizing TPS materials at these cryogenic fuel tank locations to augment 
TCAT’s TPS sizing capabilities. 

In addition, the Crank-Nicolson method was substituted for the simple implicit 
(Laasonen) Method used in the numerical solution for the interior nodal points of TCAT. 
The simple implicit method has first-order accuracy with a truncation error of 
o|^Af ^(Ax) 2 . The Crank-Nicolson Method makes use of trapezoidal differencing to 
achieve second-order accuracy with a truncation error of O^Ar^^Ax) 2 . This method 
was selected because it is the same method used in the commercial heating code SINDA. 
The Crank-Nicolson method is proven to require fewer nodes and less execution time 
than the simple impicit method. Therefore, the solution process may take less time to 
complete, and the numerical accuracy should be better than that of the simple implicit 
method. 
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1.0 Project Motivation 
Ground Hold Analysis Tool 

Conceptual design of launch vehicles is generally an iterative process in which 
disciplinary analyses are performed individually. At the Space Systems Design Lab at 
the Georgia Institute of Technology, research in the conceptual design of third generation 
launch vehicles flows through a sequence of disciplinary analyses found in a Design 
Structure Matrix. The results of each disciplinary analysis are sent to other disciplines 
whose analyses are dependent upon those results. A typical design structure matrix is 
illustrated in Figure 1 . 



Figure 1: Design Structure Matrix 

From Figure 1 we can see that the TPS sizing is dependent on the trajectory and 
configuration disciplines. More specifically, trajectory data such as altitude, relative 
velocity, and angle of attack versus time as well as geometry data are required in order to 
perform the aeroheating analysis and generate the TPS thicknesses and unit weights. 
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Once the TPS material thicknesses and unit weights are determined for the iteration, the 
data is sent to the weights and sizing discipline. 

TCAT, the Thermal Calculation Analysis Tool is a conceptual design tool 
originally developed for the Space Systems Design Lab to integrate aeroheating analysis 
into conceptual reusable launch vehicle design. TCAT is a FORTRAN code that uses 
finite difference methods to perform transient in-depth one-dimensional conduction 
analysis over the center mold line of the vehicle. TCAT couples with MINIVER, TPSA r , 
and ADS to correctly size the vehicle’s thermal protection system. MINIVER is an 
Aeroheating code that produces centerline radiation equilibrium temperatures, convective 
heating rates, and head loads over simplified vehicle geometries. TPSLT is a NASA Ames 
material properties database that is available on the World Wide Web. ADS (Automated 
Design Synthesis) is a numerical optimizer that uses algorithms to solve constrained and 
unconstrained design problems. A schematic of the TPS sizing process using TCAT is 
shown in Figure 1 . 



Figure 2: TPS Sizing Process 

When TCAT was originally created, its goal was to size TPS materials for in- 
flight aeroheating. However, sizing the TPS solely for in-flight aeroheating may not 
guarantee that the vehicle will meet all of its performance requirements. During in-flight 
operations, the surface of the vehicle experiences very high temperatures as well as high 
heat rates induced by convective heating. Therefore, TPS materials are selected so that 
the surface temperature of the materials does not exceed their operating temperature 
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limits. The materials are sized so that the temperatures experienced by the rest of the 
vehicle, namely the vehicle structure or tank structure, do not exceed some critical 
temperature, which would result in thermal failure. 

What TCAT did not account for is that during ground hold operations, the vehicle 
is held at the launch facility for a brief period while it is prepared for launch. The vehicle 
is exposed to the ambient temperature of the air, as opposed to the high temperatures 
caused by aeroheating during flight. Heat is transferred from the atmosphere to the 
surface of the TPS at a much lower rate and heat transfer occurs within the depth of the 
TPS material stack. TPS materials applied over a cryogenic tank must be sized to prevent 
the formation of ice on the outer surface of the TPS material. Since the fuel inside the 
cryogenic fuel tank is very low (-423.67 F for LH2, -297.67 for LOx), heat is extracted 
from the tank structure as well as the TPS materials directly in contact with the tank 
structure, causing temperatures to decrease throughout. If the temperature of the outer 
surface of the TPS falls below the freezing point of water (32F) then frosting will occur. 
Condensation or frosting on the surface of the TPS material is considered hazardous since 
it will negatively affect the performance of the vehicle or damage the TPS. Therefore, a 
design tool capable of sizing TPS materials during ground hold operations must be 
developed to ensure that the overall TPS thicknesses selected for the vehicle will be 
adequate to satisfy not only the in-flight aeroheating constraints, but also the ground hold 
constraints. 

Crank-Nicolson Method and TCAT 

At the conceptual design level, since the level of fidelity is not as high as it would 
be for detailed design, conceptual design tools that provide a high level of accuracy with 
low computational times are desired. TCAT is solves the one-dimensional heat equation 
numerically using a simple implicit (Laasonen) method for the interior nodes of a TPS 
material stackup. Although this method is unconditionally stable it is only first order 
accurate with a truncation error of <?|(A/)(Ar) 2 . The Crank-Nicolson method is another 
implicit numerical scheme that is also unconditionally stable; however, it makes use of 
trapezoidal differencing to achieve second-order accuracy with a truncation error of 
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Oj^Af ^.(AxJ 2 This method is the same as that used in the commercial heating code 
SINDA. The solution method of the resulting system of equations is the same as for the 
simple implicit method, but generally the Crank-Nicolson method requires fewer nodes 
and less execution time. Therefore, substitution of the Crank-Nicolson method for the 
simple implicit method may make TCAT a more valuable and practical conceptual design 
tool. 
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2.0 Ground Hold Analysis Tool Development 

The original version of TCAT operates by taking aeroheating data from 
MINIVER (a commercial aeroheating code) and sizing the TPS material for each body 
point specified by MINIVER. The data provided by MINIVER is the radiation 
equilibrium temperature and the convective heat rate versus time. During groundhold it 
is assumed that the heat transfer is constant and determined by the ambient temperature 
of the air. Since the ambient temperature of the surrounding atmosphere is constant, and 
the heat transfer rate is constant, it is not necessary to create a file like those provided by 
MINIVER. Instead, these constants can be entered explicitly into the numerical solution 
process. 

The ground hold analysis only needs to be applied to one representative body 
point above a cryogenic fuel tank location. This is because it is assumed that the 
conditions at one body point location would be the same as any other location above the 
same cryogenic fuel tank. The unit weight calculated from the single body point analysis 
can then be applied to the total area of the vehicle surface above cryogenic fuel tank 
locations for which the same TPS material will be used under the same ambient 
conditions. 

Numerical Solution Method 

TCAT uses a fully implicit method to solve the parabolic, one-dimensional 
unsteady heat conduction equation by marching in time. 


ar_ d 2 r 

dt " ° dx 2 


( 1 ) 


The boundary condition given by equation (2) is applied to the top surface of the 
TPS material. This equation is the energy balance relationship for the top surface of the 
TPS material. The convective, radiative, and conductive heat transfers are summed to 
equal zero in order to preserve the conservation of energy. 
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q^-coTS+k^-.Oax-O (2) 

However, for groundhold, a constant heat transfer mechanism is assumed, where 
the heat transfer coefficient is 1 BTU/{ft 1 ■ hour • °f). Here, the atmosphere is treated as 
an infinite atmosphere in which heat transfer to the surface TPS material can occur, but 
the temperature of the atmosphere does not change. Therefore, the boundary condition 
given in equation (3) is used for the groundhold analysis instead of equation (2). 

at x - 0 (3) 

ax 

At the backface of the TPS material the original TCAT uses the boundary 
condition in equation (4), which assumes that there is an adiabatic wall at the backface of 
the material. 


— — - 0 at jt - 1 (4) 

ax 

However, for groundhold an isothermal boundary condition is applied. 

T - T m , m const at x - L (5) 

This isothermal assumption is applied because during groundhold it is assumed 
that the temperature of the fuel remains constant, which is in direct contact with the 
backface of the cryogenic fuel tank structure. Heat can be transferred from the material 
to the fuel at the interface of the tank and fuel; however, it is assumed that the heat 
transfer that does occur does not affect the overall temperature of the fuel. This 
assumption is valid because the fuel can be treated as an infinite reservoir, due to the 
large volume of fuel that is contained inside the cryogenic fuel tanks. 
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At the surface, using the boundary condition given by equation (3) the finite 
difference equation that results is given in equation (6). 

t« . q _ 2ot|A/ t'"*i V t»*> f&\ 

• ax t; “ 4 ?"^ ■ >* T ' ' ’ 

In equation (6), the constant q is the heat transfer from the atmosphere, and is equal to the 
ambient temperature times the heat transfer coefficient. 

At the interior nodes, using a simple implicit central finite difference scheme, the 
resulting discretization is given by equation (7). 



On the back surface of the material, the isothermal boundary condition of 
equation (5) is applied and results in equation (8). 

T» - r ;* 1 - T m , - const ( 8 ) 

Examining these equations, we see that we now have a system of linear equations. 
In the numerical method used by TCAT for the aeroheating analysis, a system of non- 
linear equations resulted, because of the fourth-order temperature dependence of the 
radiative heat transfer term. Since our system is now linear, the solution process can be 
made much simpler because we no longer have to reiterate to converge to the solution. 
The system of equations can be rearranged into the form Ax - b as shown in equation (9) 
to simplify the solution process, where A is a matrix of coefficients, x is an array of 
temperatures at the next time step, and b is an array of temperatures at the current time 
step. 
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This solution of this system is easily obtained using the Thomas Algorithm, which is a 
numerical algorithm capable of solving tridiagonal banded matrices. 

Optimization Method 

For the ground hold analysis problem, the only constraint that needs to be 
enforced is the surface frosting condition. As long as the temperature of the surface of 
the TPS stays above the designated frosting condition, then the constraint is satisfied. If 
there were no heat transfer from the atmosphere, then we would expect the temperature 
profile of the surface of the TPS to decrease with time. In other words, the maximum 
temperature would be at time t » 0 and the minimum temperature would be at the end of 
the analysis. If there is heat transfer from the atmosphere, then we might expect the 
temperature of the surface to rise slightly, and then decrease to a minimum temperature at 
the end of the analysis. In both cases, we see that the minimum temperature is always at 
the end of the ground hold analysis time. Therefore, our only constraint would be to 
ensure that the temperature of the surface at the end of the analysis remains equal to or 
greater than our frosting condition temperature. Sizing the TPS to meet this constraint 
would guarantee that the design would be optimized, because it would be the minimum 
thickness required to properly insulate the surface of the TPS from the cryogenic fuel 
tank. 

The ADS optimizer used in TCAT’s aeroheating analysis is a “general problem 
optimizer,, meaning that it was created to handle many different types of engineering 
design problems. ADS has many controlling parameters that have to be fine-tuned in 


Stephen James Izon 


Pace 8 of 81 


Improvements in Thermal Protection Systems Sidne Capabilities for TCAT 

order to achieve an optimal answer. Therefore, a problem specific optimizer was written 
into the main FORTRAN code to find the optimum design for the analysis. The Newton- 
Raphson method was implemented to find the root of the function F , where F is given 
by equation (10). 


F -T* - templmt(l ) 


( 10 ) 


In this equation, T* is the surface temperature of the TPS at the final time, and 
templmt{ l) is the frosting condition temperature. Since the temperature T* is a function 
of the thickness of the material being sized, F is also a function of the same thickness. 
Therefore, to find the optimal design, we want to drive the function F to zero. 

The Newton-Raphson method is an iterative method used to find the root of a 
function. The Newton-Raphson formula is given in equation (1 1). 


x 


>>i 



01 ) 


In this notation, F is the same as in equation (10), and x is the root. In this case, we do 
not know the value of F'(^c i ), so we can approximate it using a backward finite divided 
difference. This gives us equation (12). 




( 12 ) 


Finally, substitution into equation (11) results in equation (13), which was 
integrated into the main FORTRAN code. 


*/+] 




)" F (*i) 


(13) 
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We see that two initial guesses for the root must first be made before any 
subsequent approximations for the root can be made. Therefore, in the main code, two 
guesses for the root are made before the main optimization loop begins. Once the 
difference in consecutive approximations falls below 0.0001, convergence is declared and 
the optimized thickness is obtained. 

Program Execution 

The ground hold analysis begins with the input of the ground hold time, time step 
for the analysis, ambient temperature, fuel temperature, frosting temperature, and TPS 
stackup configuration data from a web browser. Once these are input, a CGI script 
(‘single.cgi’ for a single analysis or ‘compare.cgi’ for a stackup comparison) is called. 
The CGI script takes this data and creates an input file called ‘inputs.in’ which contains 
the constraints entered by the user, and also TPS material properties based on the 
materials selected by the user and data taken from the TPSA’’ materials database. The 
FORTRAN main code (‘ground.f if the first layer is optimized or ‘ground2.f if the 
second layer is optimized) is called by the CGI script where the optimization is 
performed. The main FORTRAN code calls the subroutines ‘one.f , ‘two.f , ‘three.f , 
‘four.f , or ‘five.F, depending on the number of materials selected in the TPS material 
configuration. These subroutines solve the one-dimensional heat equation using the 
numerical method described in the previous section. Results are sent back to the main 
FORTRAN code, where the Newton-Raphson method is implemented to converge upon 
the solution. Once the solution is obtained, the results are sent back to the CGI script. 
These results are then written to the files ‘results.txt’ and ‘thickness.txt’, which are posted 
back to the user via the web interface. A schematic of the program flow logic is shown in 
Figure 3. 
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3.0 Ground Hold Analysis Test Results 

To illustrate the necessity for sizing TPS during both ground hold and in-flight 
operations, TPS materials were sized for a 10° half-angle cone. This 10° half-angle cone 
could be representative of the nosecap of a reusable launch vehicle. Figure 4 illustrates 
the geometry of the cone used to build the MINIVER input file. 



The STS-1 reentry trajectory was used as the flight profile experienced by the 
cone in the MINIVER input file. The STS-1 reentry trajectory in tabular form can be 
found in Appendix A. Body points were place two feet apart on the surface of the cone, 
resulting in 30 body points for both the leeward and windward sides. The TPS materials 
used for the analysis were CFBI, Composite Flexible Blanket Insulation, for the leeward 
surface and LI-2200, 22 lbm/ft 3 Rigid Ceramic Tile Insulation, for the windward surface. 
Schematics of the TPS material configurations are shown in Figures 5 and 6. 


Figure 5: Schematic of CFBI 
TPS Configuration 

Figure 6: Schematic of LI-2200 
TPS Configuration 


CFBI Blankets 
RTV Adhesive 
Graphite Epoxy Structure 


Ll-2200 Tiles 
RTV Adhesive 
Strain Isolator Pad 
RTV Adhesive 
Graphite Epoxy Structure 
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The results for the aeroheating analysis are given in Tables 1-3. 
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mm leeward m \ 

TPS Material 

Li-2200 Tiles 

CFBl Blankets 

TPS Unit Weight [Ib/ft 2 ] 

3.031 

0.16 


Table 3: TPS Design Results for 10* Half-Angle Cone for In-Flight Aeroheating 


The analysis was then performed using the ground hold analysis tool For each 
surface, the TPS was sized for the presence of both a liquid hydrogen (LH2) fuel tank and 
a liquid oxygen (LOx) fuel tank located inside the nosecap. The parameters of the 
analysis are given in Table 4 and the results of this analysis are given in Tables 5 and 6. 


Ground Hold Time 

7800 s (2 hours) 

i Time Step for Analysis 

10 s 

Ambient Temperature 

70 °F 

LOx Fuel Temperature 

-297.67 °F 

LH2 Fuel Temperature 

-423.67 °F 

Frosting Condition Temperature 

40 °F 


Table 4: Ground Hold Analysis Parameters 
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Table 5: TPS Sizing Results for Ground Hold Analysis Above a LOx Fuel Tank 
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3.2929 
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Table 6: TPS Sizing Results for Ground Hold Analysis Above a LH2 Fuel Tank 
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As a check to see if the results make physical sense, the temperature profiles for 
each TPS material configuration were graphed versus time. This will also give an 
indication whether or not the design constraints have been satisfied. Figures 7 and 8 
show the nodal temperature profiles for the LH2 sized TPS materials, and Figures 9 and 
10 show the nodal temperature profiles for the LOx sized TPS materials. 



Figure 7: Nodal Temperature Histories for LI-2200 Sized Over LH2 Tank 



Figure 8: Nodal Temperature Histories for CFBI Sized Over LH2 Tank 
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Figure 9: Nodal i emperature Histories for Ll-2200 Sized Over LOx Tank 
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Figure 10: Nodal Temperature Histories for CFB1 Sized Over LOx Tank 
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Improvements in Thermal Protection Systems Sizing Capabilities for TCAT 

In each of these nodal temperature histories, we see that the temperature of the 
surface never drops below 277.6 K (40 °F), which is the surface frosting constraint. 
Also, the temperature of the backface remains constant and equal to the fuel temperature 
specified for each case. Therefore, we can be confident that problems constraints have 
been met. Also, the temperature profiles are all approaching a steady state value at the 
end of the analysis. The profiles all increase in temperature through the depth of the 
stack, starting with the lowest temperature at the backface and the highest temperature at 
the surface. There is also no overlap or crossing of temperature profiles during the 
analysis. This provides assurance that the numerical method is accurately modeling the 
physical heat transfer mechanism. 

We see that for the windward surface, the unit weight of the LOx sized TPS is 
lower than the unit weight obtained from the aeroheating analysis. However, the unit 
weight of the LH2 sized TPS is higher than that of the aeroheating analysis. This means 
that the TPS sized for aeroheating is sufficiently thick to prevent surface frosting if a LOx 
tank is located in the nosecap, but it is not thick enough if a LH2 is located in the 
nosecap. 

Now, examining the results of the leeward surface, we see that the unit weights of 
both the LOx sized and LH2 sized TPS are higher than that of the aeroheating results. If 
we consider the thicknesses of the TPS alone, we can see that the TPS sized for 
aeroheating at the very first body point of the leeward side is thick enough to insulate the 
TPS from the effects of a LOx tank, but still is not thick enough for the LH2 tank. This 
means that although the TPS may be able to sustain the heat loads experienced by the in- 
flight aeroheating, the TPS still is not adequately sized for a complete mission profile. If 
the TPS results from the aeroheating analysis alone were used to construct the vehicle, 
then the majority of the surface of the vehicle would experience surface frosting. This 
would result in adverse operating conditions, which could result in TPS damage or 
possible mission failure. 

The increase in unit weights calculated from the ground hold analysis will of 
course affect the overall dry weight of the vehicle. This could possibly result in a 
significant increase in cost. Although cost is an important design driver, performance 
must not be neglected. Therefore, in order to provide the most accurate results and meet 
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all performance expectations, the TPS should be sized for both ground hold and in-flight 
aeroheating, with the TPS thicknesses and unit weights satisfying both sets of constraints 
selected for the final design. 
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4.0 Crank-Nicolson Method and TCAT 

The simple implicit method used in TCAT’s aeroheating analysis has been 
replaced with the Crank-Nicolson finite difference method in order to increase numerical 
accuracy and to decrease computational time. 

Finite Difference Equations 

The one-dimensional heat equation discretized at the interior nodes, using the 
nodal map shown in Figure 11, with a simple implict central finite difference scheme 
results in equation (14). 


,^L(r 

Ax 2 


*♦1 , 
Vi-1 ~ 2T i + 


T *it+1 Vi 

i +1 r 1 i 


(14) 



i 

Ax 

i 9 i-i 

Ax 

i • i 

f * . . , 

• n i 


Figure 11: Nodal Mesh Map for Interior Nodes 

Using the same nodal mesh map, the discretization of the one-dimensional heat 
equation with the Crank-Nicolson method results in equation (15). 



Rearranging this equation such that / is equal to the information at the next time 
step minus the information at the current time step, we arrive at equation (16). 
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This equation is now in the same form as the equations used in TCAT’s solution process. 
The resulting system of equations is iteratively solved using the Newton-Raphson 
method. First, the Jacobian Matrix is formed. 
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(17) 


The system of equations is now in the form Ax - b .given in equation (10), where A is 
the Jacobian matrix, x is the temperature change at time n + 1 and b--f at time n . 
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Solution of this system is performed by making an initial guess for the temperature at the 
next time level n + 1 and iteratively solving for the residual AT"* 1 using the Thomas 
Algorithm until convergence is achieved. 
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5.0 World Wide Web Interface 


The web version of TCAT has been completely restructured in order to provide a 
more intuitive user interface. The greeting page of TCAT welcomes the user with an 
animation in Macromedia Flash format. Following the animation, the user can select the 
‘TCAT’ link that appears in the center of the screen to go to the main content page. It is 
here that the user can select which version of TCAT to use for the analysis. 



Figure 12: Screen Shot of TCAT Welcome Page 
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Ground Hold Analysis 

If the groundhold analysis is selected, a menu appears in the content frame of the 
page where the user can choose to analyze a single TPS configuration, or compare three 
TPS configurations. If the user selects the single configuration analysis, the single 
analysis page appears in the content window. 



Figure 13: Screen Shot of TCAT Ground Hold Single Analysis Page 


The user must input the duration of the groundhold, the time step for the analysis, 
the ambient temperature of the air during groundhold, the fuel temperature for the 
isothermal boundary condition, and the surface frosting condition temperature. The user 
must also build the TPS configuration using the pull-down menus on the right side of the 
content window. While building the TPS configuration, the user can select the number of 
materials to use in the configuration. There is also the option to optimize either the first 
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or the second layer of the TPS configuration. This option is useful because some TPS 
configurations, such as RCC, SiC, or TUF1 are comprised of a thin composite surface 
layer and a diffusion layer. The surface layer remains the same thickness; however, the 
diffusion layer is sized for the analysis. 

If the user selects the multiple configuration comparison, the multiple analysis 
page appears in the content window. 



Figure 14: Screen Shot of TCAT Ground Hold Multiple Analysis Page 


Here, the user can build up to three different TPS material configurations. This 
option works in the same way as the single stackup option, the only difference being that 
there are now three sets of pull down menus in which the user can build their own 
material stackups. The results for each analysis are stored in different files, which are 
compared to determine which material stackup option provides the minimum material 
unit weight for the analysis. The results are then returned to the user in a new window. 
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Aeroheating Analysis 

The web pages for the original TCAT aeroheating program have been integrated 
into the new TCAT website. This means that both the aeroheating analysis and ground 
hold analysis can be performed at the same location. The original functionality of 
TCAT’s aeroheating analysis has remained unchanged. A screen shot of TCAT’s 
aeroheating analysis configuration is shown in Figure 15. 



MINIVER File Upload 

There is now an option available for the user to upload a MINIVER file from their 
PC to the Space Systems Design Lab’s web server. This option was made available 
because in order to perform aeroheating analyses using TCAT, a MINIVER file is 
required. To access the file upload interface, the user can simply click on the “File 
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Upload,, link in the main menu, or go to the aeroheating analysis page where a link to the 
file upload page is also available.. When the file upload screen appears, the user can then 
click on the browse button to locate a MINIVER file from their PC’s hard drive. The 
user can specify up to three files to upload at once, with a maximum file size of 1 MB. 
Once the upload button is pressed, the file is automatically uploaded to the SSDL web 
server where the aeroheating analysis files are located. The CGI script to perform the file 
upload was taken from the “Perl Services,, website, which is located at 
httD://www,DcrIserviccs.nct/en/programs/psunload/index.shtml . The html files associated 
with the ‘upload.cgi’ script have been modified to conform to TCAT’s site layout. 



Figure 16: Screen Shot of TCAT MINIVER File Upload Page 
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6.0 Conclusions 

The capabilities of TCAT at its inception, although accurate at calculating TPS 
material thicknesses when compared to commercial heating codes such as SINDA, did 
not provide complete results that are required for thermal protection system sizing 
analyses. TCAT was capable of sizing TPS materials for in-flight aeroheating; however, 
the capability to size TPS for ground hold operations was still missing. In order to meet 
this goal, TCAT’s capabilities were augmented with the addition of a ground hold 
analysis tool. This ground hold analysis tool can size the TPS materials for body points 
located above cryogenic fuel tanks using an isothermal boundary condition at the 
backface of the fuel tank structure and a constant heat transfer rate from the atmosphere 
to the surface of the TPS. The results of this analysis tool can be compared to the results 
given by the aeroheating analysis to determine the TPS thicknesses required at each body 
point location to meet both the aeroheating and ground hold temperature constraints. 

In our test case, we discovered that in some instances the TPS sized for 
aeroheating may be sufficiently thick to prevent surface frosting. This was evidenced in 
the case of the TPS sized for the windward side of the 10° half-angle cone with a LOx 
tank located inside the cone. However, for the rest of the analyses it was determined that 
the TPS sized for aeroheating were not sufficient to meet the constraints of the ground 
hold. Therefore, the design tool created specifically for the sizing of TPS during ground 
hold should be integrated into the TPS design process. 

Even with the ground hold analysis tool developed, TCAT’s capabilities are not 
yet entirely developed. A method of integrating the ground hold and aeroheating analysis 
should be developed so that the TPS thicknesses of each individual analysis do not need 
to be compared on a point-by-point basis. This could be achieved by integrating the 
ground hold temperature and heat transfer profiles for each body point into the MINIVER 
file that includes the radiation equilibrium temperature and heat transfer profiles during 
in-flight trajectory. The entire analysis could be performed using this single input file, 
with the method of analysis switching from ground hold constraints to in-flight 
aeroheating constraints, returning the optimal TPS thickness that satisfies both sets of 
constraints to the user. 
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Appendix A: STS-1 Reentry Trajectory 


Time (s) 

Attitude (ft) 

Velocity (ft/s) 

Angle of Attack 

Bank Angle 

45.3 

373800 

24590 

41.26 

0 

90.3 

351500 

24620 

41.21 

0 

135.3 

329500 

24640 

40.46 

0 

180.3 

308000 

24660 

40.65 

0 

225.3 

288000 

24660 

41.9 

0 

270.3 

269700 

24610 

39.39 

0 

315.3 

256300 

24470 

40.71 

0 

360.3 

250200 

24220 

41.74 

0 

405.3 

247000 

23920 

39.95 

0 

450.3 

244600 

23610 

39.28 

0 

495.3 

242500 

23280 

39.66 

0 

545.3 

240000 

22880 

39.07 

0 

569.3 

238800 

22680 

39.32 

0 

593.3 

237700 

22470 

39.02 

0 

617.3 

236600 

22250 

39.25 

0 

641.3 

235300 

22020 

39.48 

0 

665.3 

233700 

21780 

39.97 

0 

689.3 

232000 

21530 

40.01 

0 

713.3 

232300 

21260 

40.25 

0 

737.3 

230300 

20980 

40.46 

0 

761.3 

227900 

20680 

40.14 

0 

785.3 

225200 

20360 

40.16 

0 

809.3 

223000 

20020 

40.16 

0 

833.3 

219000 

19650 

40.27 

0 

857.3 

215300 

19240 

40.4 

0 

881.3 

211400 

18790 

40.12 

0 

905.3 

207600 

18300 

40.3 

0 

929.3 

205600 

17760 

40.01 

0 

953.3 

202500 

17180 

42 

0 

977.3 

197200 

16540 

40.9 

0 

1001 

192200 

15870 

40.75 

0 

1025 

187200 

15130 

39.87 

0 

1049 

182600 

14350 

39.43 

0 

1073 

179600 

13570 

39.53 

0 

1097 

176600 

12790 

39.75 

0 

1121 

172500 

12010 

38.96 

0 

1145 

167500 

11220 

38.12 

0 

1169 

161900 

10440 

36.97 

0 

1240 

150000 

8336 

34.07 

0 

1302 

133900 

6757 

28.2 

0 

1364 

117400 

5342 

23.05 

0 

1426 

106200 

4064 

20.27 

0 

1488 

89100 

2917 

16.47 

0 

1550 

76350 

1915 

10.86 

0 

1612 

57440 

1151 

7.79 

0 

1674 

41670 

800 

7.72 

0 

1736 

27600 

681 

7.09 

0 

1860 

3337 

512 

3.78 

0 

1925 

-3 

202 

-12.2 

0 
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Appendix B: Windward Aerobeating Analysis Input File 


10 

48 

12 

12 

12 

12 

0.899 

0.074420 

0.396700 

0.042000 

0.396700 

0.601300 

352.400000 

1410.000000 

86.500600 

1410.000000 

1576.232900 

0.628000 

1.168000 

1.318000 

1.168000 

0.833300 

294.2611111 

0.075000 

0.003000 

0.004000 

0.003000 

0.003000 

277.5944444 

1644 

561.100000 

561.100000 

505.222000 

505.222000 

561.100000 

561.100000 

588.889 

19.99999999 

6 

5 

5 

1340 


Stephen James 


! /time step 

! /number of nodes of material 1 
! /number of nodes of material 2 
! /number of nodes of material 3 
! /number of nodes of material 4 
! /number of nodes of material 5 
!/emissivity of material 1 
! /thermal conductivity of material 1 

! /thermal conductivity of material 2 

! /thermal conductivity of material 3 

! /thermal conductivity of material 4 

! /thermal conductivity of material 5 

! /density of material 1 
! /density of material 2 
! /density of material 3 
! /density of material 4 
! /density of material 5 
! /specific heat of material 1 
! /specific heat of material 2 
! /specific heat of material 3 
! /specific heat of material 4 
'/specific heat of material 5 
! /initial temp (deg K) 

! /TPS thickness of material 1 
! /TPS thickness of material 2 
! /TPS thickness of material 3 
! /TPS thickness of material 4 
! /TPS thickness of material 5 
! /Temperature limit of material 1 
! /Temperature limit of material 1 
! /Temperature limit of material 2 
! /Temperature limit of material 2 
! /Temperature limit of material 3 
! /Temperature limit of material 3 
! /Temperature limit of material 4 
! /Temperature limit of material 4 
! /Temperature limit of material 5 
! /backf ace temperature limit 
!/istrat parameter 
!/iopt parameter 
!/ioned parameter 
!/iprint parameter 
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Appendix C: Leeward Aeroheating Analysis Input File 


3 

! /number of materials 


10. 

!/time step 


40 

! /number of nodes in material 1 


10 

'/number of nodes in material 2 


10 

’/number of nodes in material 3 


0.8856 

’/emissivity of material 1 


0.03821 

’/thermal conductivity of material 

1 

0.3967 

’/thermal conductivity of material 

2 

0.6013 

! /thermal conductivity of material 

3 

96.118 

! /density of material 1 


1410 

'/density of material 2 


1576.2329 

! /density of material 3 


0.7821 

! /specific heat of material 1 


1.168 

! /specific heat of material 2 


0.8333 

! /specific heat of material 3 


300. 

’/initial temp (deg K) 


0.05 

! /TPS thickness of material 1 (meters) 

0.003 

’/TPS thickness of material 2 (meters) 

0.003 

! /TPS thickness of material 3 (meters) 

1477.44 

’/Temperature limit of material 1 


1477.44 

! /Temperature limit of material 1 


561.1 

! /Temperature limit of material 2 


561.1 

! /Temperature limit of material 2 


588.889 

! /Temperature limit of material 3 


421.889 

! /backf ace temperature limit 


6 

!/istrat parameter 


5 

!/iopt parameter 


5 

!/ioned parameter 


0000 

! /iprint parameter 
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Appendix D: 

5 

10 

48 

12 

12 

12 

12 

0.899 

0.074420 

0.396700 

0.042000 

0.396700 

0.601300 

352.400000 

1410.000000 

86.500600 

1410.000000 

1576.232900 

0.628000 

1.168000 

1.318000 

1.168000 

0.833300 

294.2611111 

0.075000 

0.003000 

0.004000 

0.003000 

0.003000 

277.5944444 

1644 

561.100000 

561.100000 

505.222000 

505.222000 

561.100000 

561.100000 

588.889 

90 

6 
5 
5 

1340 


TQA1 


Windward Ground Hold Analysis for LOx Tank Input File 

!/ number of materials 
!/time step 

! /number of nodes of material 1 
! /number of nodes of material 2 
! /number of nodes of material 3 
! /number of nodes of material 4 
! /number of nodes of material 5 
!/emissivity of material 1 
! /thermal conductivity of material 1 
.'/thermal conductivity of material 2 
! /thermal conductivity of material 3 
! /thermal conductivity of material 4 
! /thermal conductivity of material 5 
! /density of material 1 
! /density of material 2 
! /density of material 3 
! /density of material 4 
! /density of material 5 
! /specific heat of material 1 
! /specific heat of material 2 
!/ specific heat of material 3 
! /specific heat of material 4 
! /specific heat of material 5 
! /initial temp (deg K) 

!/TPS thickness of material 1 
\ /TPS thickness of material 2 
! /TPS thickness of material 3 
! /TPS thickness of material 4 
! /TPS thickness of material 5 
! /Temperature limit of material 1 

.'/Temperature limit of material 1 

! /Temperature limit of material 2 

! /Temperature limit of material 2 

! /Temperature limit of material 3 

! /Temperature limit of material 3 

! /Temperature limit of material 4 

! /Temperature limit of material 4 

! /Temperature limit of material 5 

! /backf ace temperature limit 
!/istrat parameter 
!/iopt parameter 
!/ioned parameter 
! /iprint parameter 
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3 

10 

40 

10 

10 

0.8856 

0.038210 

0.396700 

0.601300 

96.118000 

1410.000000 

1576.232900 

0.782100 

1.168000 

0.833300 

294.2611111 

0.050000 

0.003000 

0.003000 

277.5944444 

1477.44 

561.100000 

561.100000 

588.889 

90 

6 

5 

5 

1340 
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Leeward Ground Hold Analysis for LOx Tank Input File 

! /number of materials 
!/time step 

! /number of nodes of material 1 
! /number of nodes of material 2 
! /number of nodes of material 3 
! /emissivity of material 1 
! /thermal conductivity of material 1 

! /thermal conductivity of material 2 

! /thermal conductivity of material 3 

! /density of material 1 
! /density of material 2 
! /density of material 3 
’/specific heat of material 1 
! /specific heat of material 2 
! /specific heat of material 3 
! /initial temp {deg K) 

! /TPS thickness of material 1 
! /TPS thickness of material 2 
! /TPS thickness of material 3 
’/Temperature limit of material 1 

! /Temperature limit of material 1 

! /Temperature limit of material 2 

! /Temperature limit of material 2 

! /Temperature limit of material 3 

!/backface temperature limit 
!/istrat parameter 
! /iopt parameter 
!/ioned parameter 
’/iprint parameter 
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Appendix F: 

5 

10 

48 

12 

12 

12 

12 

0.899 

0.074420 

0.396700 

0.042000 

0.396700 

0.601300 

352.400000 

1410.000000 

86.500600 

1410.000000 

1576.232900 

0.628000 

1.168000 

1.318000 

1.168000 

0.833300 

294.2611111 

0.075000 

0.003000 

0.004000 

0.003000 

0.003000 

277.5944444 

1644 

561.100000 

561.100000 

505.222000 

505.222000 

561.100000 

561.100000 

588.889 

19.99999999 

6 
5 
5 

1340 
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Windward Ground Hold Analysis for LH2 Tank Input File 

! /number of materials 
! /time step 

! /number of nodes of material 1 
! /number of nodes of material 2 
! /number of nodes of material 3 
! /number of nodes of material 4 
) /number of nodes of material 5 
!/emissivity of material 1 
! /thermal conductivity of material 1 

! /thermal conductivity of material 2 

! /thermal conductivity of material 3 

! /thermal conductivity of material 4 

! /thermal conductivity of material 5 

! /density of material 1 
! /density of material 2 
! /density of material 3 
! /density of material 4 
! /density of material 5 
1 /specific heat of material 1 
! /specific heat of material 2 
! /specific heat of material 3 
! /specific heat of material 4 
! /specific heat of material 5 
! /initial temp (deg K) 

! /TPS thickness of material 1 
! /TPS thickness of material 2 
! /TPS thickness of material 3 
! /TPS thickness of material 4 
! /TPS thickness of material 5 
! /Temperature limit of material 1 

! /Temperature limit of material 1 

! /Temperature limit of material 2 

! /Temperature limit of material 2 

[/Temperature limit of material 3 

! /Temperature limit of material 3 

[/Temperature limit of material 4 

I/Temperature limit of material 4 

! /Temperature limit of material 5 

I /backface temperature limit 
I /istrat parameter 
! /iopt parameter 
! /ioned parameter 
! /iprint parameter 
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Appendix G: Leeward Ground Hold Analysis for LH2 Tank Input File 


3 

10 

40 

10 

10 

0.8856 

0.038210 

0.396700 

0.601300 

96.118000 

1410.000000 

1576.232900 

0.782100 

1.168000 

0.833300 

294.2611111 

0.050000 

0.003000 

0.003000 

277.5944444 

1477.44 

561.100000 

561.100000 

588.889 

19.99999999 

6 

5 

5 

1340 


/number of materials 
/time step 

/number of nodes of material 1 
/number of nodes of material 2 
/number of nodes of material 3 
/emissivity of material 1 
/thermal conductivity of material 1 

/thermal conductivity of material 2 

/thermal conductivity of material 3 

/density of material 1 
/density of material 2 
/density of material 3 
/specific heat of material 1 
/specific heat of material 2 
/specific heat of material 3 
/initial temp (deg K) 

/TPS thickness of material 1 
/TPS thickness of material 2 
/TPS thickness of material 3 
/Temperature limit of material 1 
/Temperature limit of material 1 
/Temperature limit of material 2 
/Temperature limit of material 2 
/Temperature limit of material 3 
/backface temperature limit 
/istrat parameter 
/iopt parameter 
/ioned parameter 
/iprint parameter 
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Appendix H: ‘single.cgi’ Ground Hold Single Analysis CGI Script 


# ! /usr/local/bin/perl 

read (STDIN, Sbuffer, $ENV { * CONTENT_LENGTH ' > ) ; 

{•pairs * split (/t/, $buffer); 
foreach $pair (gpairs) 

{ 

<$name, $value) - split (/-/, $pair) ; 

$value — tr/+/ /; 

$value — s/I { [a-fA-FO-9] (a-fA-FO-9) ) /pack("C", hex($l))/eg; 
$value — • s/- ! / ~!/g: 

$FORM($name) - $ value; 

) 

* — - - 

# assign information from input website to appropriate variable 

# — 


$matnum«$FORM(matnum) ; 
$method-$FORM (method) ; 
$material_l«$FORM(material_l ) 
$material~2«$FORM(material~2 ) 
$material_3-$FORM(material_3 ) 
$material_4«$FORM(material_4 ) 
$material_5-$FORM(material_5 ) 
$holdtime-$ FORM {holdtime } ; 

$ t ime step-$ FORM { times tep) ; 
$ambient-$FORM { ambient ) ; 

$f uel temp- $ FORM { fuel temp } ; 

$ frost- $ FORM (frost) ; 


* — 

# obtain inputs from command line 

# 

# List of variables 

# $bodyarea - surface area of the body in square feet 

# Sbodytile - fraction of bodyarea covered with tiles 

# $bodyblanket - fraction of bodyarea covered with blankets 
I $filename - name of MINIVER 1_<>, s file 

# Ssharp - indicates the use of SHARP materials (yes or no) 

I Sbodyside - vehicle side analysis is conducted on (leeward or windward) 

# $material - type of material to be used along with backface material 

# — - 

#**#*#*##*•###•• t ########!###♦###•# I##### ###########♦####♦ !##*#*### I ######*#♦###### 
##### 

# Prompt the user for inputs. 

# 

##t######t##M####f#f#########f#############t#l###t#######tt#t########Mt##MM###M#f### 

««**« 

l$holdtime-7800; 

#$timestep«10; 

♦$heatrate-l; 

♦$ambient-70; 

#$fueltemp— 423.67; 

♦ $ frost-40; 

Iprint "Enter the amount of time, in seconds, the vehicle will be in groundhold . \n" ; 
#$holdtime - readline(*STDIN); 

Ichop ($holdtime) ; 

Iprint "Enter the timestep, in seconds, for the analysis . \n"; 

#$timestep - readline (*STDIN) ; 

Ichop (Stimestep) ; 

♦print "Enter the heat transfer coefficient (BTU/ (f t A 2*hour*F) ) . \n" ; 

#$heatrate - readl ine ( *STDIN) ; 

♦ chop ($heatrate) ; 
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♦print "Enter the ambient temperature (F) of the groundhold location . \n" ; 

#$ambient - readline (*STDIN) ; 

♦chop ($ambient) ; 

♦print "Enter the temperature [F] of the cryogenic fuel (ie. the isothermal boundary 
conidtion) . \n"; 

♦$fueltemp - readline (*STDIN) ; 

♦chop {Sfueltemp) ; 

♦print "Enter the frosting condition Temperature [F] of the surface TPS material . \n"; 
♦$frost ■ readline (*STDIN) ; 

♦chop ($frost) ; 

open (TIME, ">time") ; 

print TIME "$holdtime ! /Time of Groundhold"; 

close (TIME); 

$heatrate - 1; 

$heattransfer - $heatrate • $ambient *3.15459074506; 
open (HEATXFER, ">heat") ; 

print HEATXFER "$heattransf er !/Heat Flux"; 

close (HEATXFER) ; 

Sambient - ($ambient-32 ) * (5/9) +273 . 15; 

$f rost - ($frost-32)* (5/9) +273. 15; 

Sfueltemp - ($fueltemp-32) * (5/9) +273. 15; 
open (TANK, ">fueltemp") ; 

print TANK "$fueltemp ! /Temperature of fuel"; 

close (TANK) ; 

♦print "Please enterXn"; 

♦print " 1. if you would like to create a new TPS material stackup for the 

analysis, or\n"; 

♦print " 2. if you would like to use an existing TPS material stackupNn"; 

♦$material_option - readl ine ( *STDIN) ; 

♦chop ($material_option) ; 

♦print "Please enter the number of materials in the TPS stackup (up to 5).\n"; 
♦Smatnum - readline ( *STDIN) ; 

♦chop (Smatnum) ; 

♦if ($material_option eq 1 ) ( 

♦ for <$loopnum«l; Sloopnum <- Smatnum; $loopnum++) ( 

♦ if (Sloopnum eq 1) { 

♦ print "Please select the surface material of the TPS stackup. \n"; 

♦ print "\n"; 

♦ print " Laminate Surfaces \n"; 

♦ print "RCC_Laminate_Surf ace\n"; 

♦ print "SiC_Laminate_Surface\n"; 

♦ print "TUFI_Composite_Surface\n"; 

♦ print " Tiles . \n" ; 

♦ print "TUFI_AETB8_Dif fusion_Layer\n"; 

♦ print "TUFI_AETB12_Dif fusion_Layer\n"; 

♦ print "FRCI 12\n" ; 

♦ print "FRCl20\n"; 

♦ print "LI900\n"; 

♦ print "Ll2200\n" ; 

♦ print " Blankets \n"; 

♦ print "AFRSI\n"; 

♦ print "CFBI\n"; 

♦ print "DURAFRSIXn"; 

♦ print "PBlNn"; 
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print "TABlXn"; 
print "Xn"; 

print "Enter the material name as it appears above. Xn"; 
$material_l * readline ( *STDIN) ; 
chop ($material_l ) ; 

) elsif ($loopnum eq 2) { 

print "Please select second material of the TPS stackupAn"; 
print "\n"; 

print " Tiles An"; 

print "TUFI_AETB8_Diffusion_Layer\n"; 

print "TUFI_AETB12 Dif f usion_Layer\n" ? 

print "FRCI12\n"; 

print "FRCI20\n"; 

print "Ll900\n"; 

print "Ll2200\n"; 

print "RCC\n"; 

print "SiC\n"; 

print " Blankets \n"; 

print "AFRSIXn"; 
print "CFBIXn"; 
print "DURAFRSlXn"; 
print "PBI\n"; 
print "TABlXn"; 

print " Intermediates Xn"; 

print "Strain_Isolator_Pad\n"; 
print "RTV_Adhesive\n" ; 
print "Rohacell_Foam\n"; 

print " Structures Xn"; 

print "GrEx_Structure\n"; 
print "Al_Structure\n"; 
print "TiAl_Structure\n"; 
print "Xn"; 

print "Enter the material name as it appears above An"; 
$material_2 - readline (*STDIN) ; 
chop ($ma ter ial_2) ; 

} elsif (Sloopnum eq 3> { 

print "Please select third material of the TPS stackupAn"; 
print "Xn"; 

print " Tiles An"; 

print "TUFI_AETB8_Diffusion_Layer\n"; 

print "TUFI_AETB12 Dif f usion_Layer\n" ; 

print "FRCI12Xn"; 

print "FRCI20\n" ; 

print "Ll900\n"; 

print "Ll2200\n" ; 

print "RCC\n"; 

print "SiC\n"; 

print " Blankets \n"; 

print "AFRSIXn"; 
print "CFBIXn"; 
print "DURAFRSlXn"; 
print "PBlXn"; 
print "TABlXn"; 

print " Intermediates \n"; 

print "Strain Isolator_Pad\n"; 
print "RTVJVdhesiveXn"; 
print "Rohacell_Foam\n"; 

print " Structures \n"; 

print "GrEx_StructureXn"; 
print "Al_Structure\n"; 
print "TiAl_StructureXn"; 
print "\n"; 

print "Enter the material name as it appears above An"; 
$mater ial_3 - readl ine ( *STDIN) ; 
chop ($material_3) ; 

# ) elsif (Sloopnum eq 4) { 
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print "Please select fourth material of the TPS stackup An"; 
print " \n"; 

print " Tiles An"; 

print "TUFI_AETB8_Diffusion_Layer\n"; 

print "TUFI AETB12_Diff usion_Layer\n"; 

print "FRCI12\n"; 

print "FRCI20\n"; 

print "LI900\n"; 

print "LI2200\n"; 

print "RCC\n"; 

print "SiC\n"; 

print " Blankets \n"; 

print "AFRSINn"? 
print "CFBI\n"? 
print "DURAFRSIXn"; 
print "PBI\n"; 
print *TABI\n"; 

print " Intermediates \n"; 

print "Strain_Isolator_Pad\n"; 
print "RTV_Adhesive\n"; 
print "Rohacell_Foam\n"; 

print " Structures \n"; 

print "GrEx_Structure\n"; 
print "Al_Structure\n"; 
print "TiAl_Structure\n"; 
print "\n"; 

print "Enter the material name as it appears above An"; 
$material_4 - readline (*STDIN) ; 
chop ($material_4 ) ; 


) elsif <$loopnum eq 5) ( 


print "Please select fifth material of the TPS stackup An"; 
print "\n"; 

print " Tiles An"; 

print "TUFI_AETB8_Dif fusion_Layer\n"; 

print "TUFI_AETB12_Diff usion_Layer\n"; 

print "FRCI12\n"; 

print "FRCI20\n"; 

print "LI 900\n"; 

print "LI2200\n"; 

print "RCC\n"; 

print "SiC\n"; 

print " Blankets \n"; 

print "AFRSI\n"; 
print "CFBI \n" ; 
print "DURAFRSI\n"; 
print "PBI\n"; 
print "TABI\n"; 

print " Intermediates \n"; 

print "Strain_Isolator_Pad\n"; 
print "RTV_Adhesive\n"; 
print "Rohacell_Foam\n"; 

print " Structures \n"; 

print "GrEx_Structure\n"; 
print "Al_Structure\n"; 
print "TiAl_Structure\n"; 
print "\n"; 

print "Enter the material name as it appears above An"; 
$material_5 - readline (*STDIN) ; 
chop <$material_5) ; 


) 


) 


) 

print "\n"; 

print "Please enter the optimization methodAn"; . 
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♦print " 1. If you would like to optimize the surface material thickness, or\n"; 
♦print " 2. If you would like to optimize the second TPS material thickness . \n" ; 
♦print "\n"; 

♦print "NOTE: If either of the laminate materials were selected as the\n"; 

♦print " surface TPS material, the second optimization method will\n"; 

♦print " be automatically selected . \n"; 

♦$method - readline (*STDIN) ; 

♦chop (Smethod) ; 

♦print "\n"; 

$material_option-l ; 

♦$matnum»5; 

♦Smaterial_l«AFRSI2500; 

♦Smaterial”l-TUFI_A£TB12_Diffusion_Layer; 

♦ $ mater ial__2-RTV_Adhesive; 

♦Smaterial_3-GrEx_Structure; 

#$material_3-Strain_Isolator_Pad; 

♦ $mat er i a 1_4 -Rohacel l_Foam; 

♦ $material_5-GrEx__Structure; 

if ($material_option eq 1) { 

for ($loopnum*l; Sloopnum <» Smatnum; $loopnum++) { 
if ($loopnum eq 1) { 

open (DAT, "<$material_l") ; 

$i - 0; 

while (<DAT>) { 

Si - Si + 1; 
my (5 line) - 
chomp (Sline) ; 

if ($i eq 1) { 

Snodes ( 1 ] - Sline; 

} elsif (Si eq 2) { 

Scond [ 1 J * Sline; 

) elsif (Si eq 3) < 

Sdensityll] - Sline; 

) elsif (Si eq A) { 

$specheat(l] - Sline; 

) elsif (Si eq 5) ( 

Sthick { X ) - Sline; 

) elsif (Si eq 6) ( 

$ tempi imit (1 1 -Sline; 

> elsif (Si eq 7) ( 

Semissivity-Sline; 

} 

> 

close (DAT) ; 

) elsif (Sloopnum eq 2) ( 

open (DAT, "<$mater ial__2") ; 

Si - 0; 

while (<DAT>) { 

Si - Si + 1; 
my (Sline) - $_? 
chomp (Sline) ; 

if (Si eq 1) { 

Snodes (2 J - Sline; 

} elsif (Si eq 2) { 

Scond ! 2 ] - Sline; 

) elsif (Si eq 3) { 
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$density[2) * Sline; 

) elsif ($i eq 4) ( 
$specheat[2] - Sline; 
} clsif ($i eq 5) ( 
$thick[2] - Sline; 

} elsif (Si eq 6) { 
Stemplimit {2 J«$line; 

> 


close (DAT) ; 

} elsif (Sloopnum eq 3) { 
open (DAT, "<$material_3") ; 

$i - 0; 

while (<DAT» ( 

Si - $i ♦ 1? 
my (Sline) - $_; 
chomp (Sline) ; 

if (Si eq 1) ( 

$ nodes (3) • Sline; 

} elsif (Si eq 2) { 
$cond(3) * Sline; 

) elsif (Si eq 3) ( 
Sdensity[3] - Sline; 

) elsif (Si eq 4) ( 
$specheat(3] - Sline; 
) elsif (Si eq 5) ( 
Sthick (3) - Sline; 

) elsif (Si eq 6) { 

S tempi imit 1 3] -Sline; 

) 


close (DAT) ; 

} elsif (Sloopnum eq 4) ( 
open (DAT, r, <$material__4 w ) ; 

Si - 0; 

while «DAT>) l 
Si - Si + 1; 
my (Sline) - $_; 
chomp (Sline) ; 

if (Si eq 1) { 

$nodes[4J - Sline; 

) elsif (Si eq 2) { 
$cond[4] - Sline; 

) elsif (Si eq 3) ( 
$density(4] - Sline; 

> elsif (Si eq 4) ( 
$specheat(4) - Sline; 
) elsif (Si eq 5) { 
Sthick (4] - Sline; 

) elsif (Si eq 6) { 

S tempi imit (4] -Sline; 

) 


close (DAT) ; 
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I elsif ($loopnum eq 5) { 
open (DAT, "<$material_5") ; 

$i - 0; 

while (<DAT>) { 

$i - $i + 1 ; 
my <$line) - 
chomp ($line) ; 

if (Si eq 1) { 

$nodes[5) - $line; 

) elsif ($i eq 2) { 
$cond [5] - $line; 

} elsif ($i eq 3) { 
$densityl5] - $line; 

) elsif ($i eq 4) { 

$specheat[5] - $line; 
) elsif ($i eq 5) { 
$thick [5] - $line; 

} elsif ($i eq 6) [ 
$templimit [ 5) -$line; 

} 


close (DAT) ; 


) 


) 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦#♦#♦♦#♦#♦♦ 

»##*#*#*#» 

# Create the TPS stackup file by reading the properties of each material selected and 
entering # 

# their values into a input file recognizable by T-CAT. 

# 

ft**###### 


♦ system ("touch stackup. in" ) ; 

♦ system ("chmod 777 stackup . in" ) ; 
open (STACKUP, ">stackup . in" ) ; 
print STACKUP "Smatnum 

print STACKUP "$timestep 


♦create TPS stackup file 

♦change permissions of TPS stackup file 

! /number of materials\n"; 

! /time step\n"; 


for ($loopnum-l; Sloopnum <■ Smatnum; $loopnum++) { 

printf STACKUP "li ! /number of nodes of material 

$loopnum\n", $nodes [$loopnum] ? 

) 

print STACKUP "$emissivity ! /emissivity of material l\n"; 

for ($loopnum-l ; $loopnum <- Smatnum; $loopnum++) { 

printf STACKUP "%f ! /thermal conductivity of material 

$loopnum\n", $cond { Sloopnum) ; 

I 

for ($loopnum*l; $loopnum <« Smatnum; $loopnum++) ( 

printf STACKUP "%f ! /density of material $loopnum\n", $density ( $loopnum) ; 

> 

for ($loopnum«l; $loopnum <- Smatnum; $loopnum++) ( 

printf STACKUP "%f ! /specific heat of material 

$loopnum\n", Sspecheat I $loopnum] ; 

} 

print STACKUP "$ambient ! /initial temp (deg K)\n"; 
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for ($loopnum«l; $loopnum <« $matnum; $loopnum++) { 

printf STACKUP "%f ! /TPS thickness of material $loopnum\n" , $ thick I $loopnum] ; 

> 


print STACKUP "$frost ! /Temperature limit of material l\n"; 

print STACKUP "$templimit ( 1 J ’/Temperature limit of material l\n"; 

for ($loopnum«2; Sloopnum <* $matnum-l; $loopnum++) { 

printf STACKUP *%f ! /Temperature limit of material 

$loopnum\n", $templimit ($loopnum] ; 

printf STACKUP "%f ’/Temperature limit of material 

$loopnum\n", $templimit l$loopnum] ; 

) 


print STACKUP 
print STACKUP 
print STACKUP 
print STACKUP 
print STACKUP 
print STACKUP 


*$templimit [ $matnum) 
'$fueltemp 
•6 
*5 
*5 

*1340 


’/Temperature limit of material $matnum\n"; 
! /backface temperature limit\n**; 

’/istrat parameter\n"; 

! / iopt parameter \n"; 

!/ioned parameter \n"; 

! /iprint parameter \n"; 


close (STACKUP) ; 

# system ("chmod 777 stackup . in" ) ; 


) else ( 


print "Enter the material for TPS An"; 

print " Tiles An"; 

print "AETB8_Al\n"; 
print "AETB8 GrEx\n"; 
print "AETB12 Al\n"; 
print "AETB12~GrEx\n" ; 
print "FRCI12_Al\n"; 
print "FRCI12 GrEx\n"; 
print "FRCI20~Al\n"; 
print "FRCI20_GrEx\n"; 
print "Ll900_Al\n"; 
print "LI 900 GrEx\n"; 
print "Ll2200_Al\n"; 
print "Ll2200_GrEx\n"; 
print "RCC_Al\n"; 
print "RCC_GrEx\n"; 
print "SiC_Al\n"; 
print "SiC_GrEx\n"; 
print "RCC_Al\n" ; 
print "RCC_GrEx\n" ; 
print "TUFI_Al\n"; 
print "TUFI_GrEx\n"; 

print " Blankets \n"; 

print "AFRSI_Al\n"; 
print "AFRSI GrEx\n"; 
print "CFBI Al\n"; 
print "CFBI~GrEx\n"; 
print "DURAFRSI_Al\n"; 
print "DURAFRSI~GrEx\n"; 
print "PBI_Al\n"; 
print "PBI_GrEx\n"; 
print "TABI_Al\n"; 
print "TABI_GrEx\n"; 
print "\n"? 

print "Enter the material name as it appears above An"; 
$material * readline ( *STDIN) ; 
chop ($material) ; 


1 

########### 
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44444I4444444444444444444944444444II44444444444444444444444444444444444444444444444444444 

444444444444 

$numtimes * $holdtime/$timestep ♦ 1; 

Stime [1] - 0; 

for ($i -2; $i <- $numtimes; $i++) { 

$time{$i] - $time[$i-l] 4 $timestep; 


) 

♦system {"touch l_stackup.s") ; 4create mock MINIVER file 

♦system ("chmod 777 l_stackup. s") ; 4change permissions of file 
open (MINI, *>l_stackup. s"j ; 

print MINI " Mock MINIVER file for input into ground hold analysis 10 

1 . 0\n"; 

for ($i«l; $i <- $numtimes; $i+4) { 

print MINI "$time[$i) $ambient $heattransfer\n"; 


) 

print MINI "-100.0"; 
close (MINI) ; 


♦ - - 

♦ MINIVER file is partitioned into body point files 

♦ partition the MINIVER file into separate body point files that include time, 

♦ heat rate, and radiation equilibrium temperature 

♦ — 


$i«0; ^initialization of counter index "i" 

$f ilecount-0; 

open (FILE, M <l_stackup.s") ; 4open statement with filehandle FILE 


until (eof FILE) { 
reached 

START: $_-<FILE>; 
if </ A . (73) (\S* ) / ) 
$title-$l; 


♦goes through the MINVER file until the end is 

♦skips the first line of the data file 

♦searches for the bodypoint label 


4 system {"touch groundhold") ; 4create body point file 

4 system("chmod 111 groundhold"); 4change permissions inorder to write to body point 

file 

$f ilecount*$f ilecount+l; 

) 4 end of if statement 

while (<FILE>) { 4while loop that reads thru each line of the file 

Si-$i+l; 

/*\sM\S*)\s*<\S*)\s* (\S*)\s*(\S*)\s*(\S*)\s*(\S*)\s*(\S*)\s* (\S*) Vs* (\S*) \sM\S* > \s* (\S* 
) \s* (\S* ) \s* (\S*) /; 

$time($i J-$l; 4array declaration for time 

$tradeq[$i)«$2; 4array declaration for rad eq temp 

$qconv[$i)«$3; 4array declaration for conv heat rate 

$n-$i; 

if ($1—100) { 

open (NEWFILE, ">groundhold") ; 4open statement for new file "newsample" 

$ times teps»$n-l / 

printf NEWFILE "$timesteps \n"; ♦ formatted print statement" 

for ($i-l; $i<$n; $i44) { 

printf NEWFILE "If \n" , $time ( $i ) ; 
printf NEWFILE "If \n" , $qconv($i) ; 
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print f NEWFILE "If \n" , $tradeq [ $i] ; 

) fend of for loop 

close (NEWFILE); 

$i»0; 

goto START 

) fend of if statement 

) fend of while loop 
) fend of until loop 
close (FILE) ; 

f - — — 

f incrementally move each body point file into "miniver . in" and execute teat 
f perl script that creates the filenames needed inorder to batch execute teat 
f filename are created via input from the user 
f inputs (file_start number, file_end_number, body_side) 

f — - 

$start-l; fstart of the file index 

f Send* ( S f ilecount-1) /2? fend of the file index 

$end-l; 

open (FILE, ">inputs_for outputs") ; 
printf FILE "l\n"; 
close (FILE); 

f system("touch thickness . txt") ; 
f system ( "chmod 777 thickness.txt"); 
open (FINALOUTPUT,">thickness. txt") ; 

print FINALOUTPUT "TPS Ground hold design results\n\n"; 

f — 

f system ("touch miniver . in") ; 

f system ("chmod 777 miniver . in") ; 

system ("cp groundhold miniver . in") ; 

I 

f assign material properties 

I — 

if ($material_option eq 1) ( 

Sinputfile * "stackup. in"; 

) else ( 

if (Smaterial eq "AETB8_A1") ( 

Sinputf ile-"AETB 8_AL 5inputs.in"; 

) 

if ($material eq "AETB8_GrEx" ) ( 

$ input file“"AETB_8_GrEx_5 inputs . in"; 

) 

if ($mater ial eq "AETB12_A1") ( 

$ input file«"AETB 12_AL_5inputs . in"; 

} 

if ($material eq "AETB1 2_GrEx" ) { 

$inputf ile»"AETB_12_GrEX_5inputs . in"; 

) 

if (Smaterial eq "FRCI12 Al") { 

$inputfile»"FRCI_12_AL_5inputs.i"; 

) 

if (Smaterial eq "FRCI 12_GrEx") ( 

Sinputf ile-"FRCI_12_GrEx_5inputs . in"; 

1 

if (Smaterial eq "FRCI20_A1") { 

Sinputf ile-"FRCI_20_AL_5inputs . in"; 

) 

if (Smaterial eq " FRCl20_GrEx") ( 

Sinput f ile=” FRCI_20_GrEx_5inputs . in"; 
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if (Smaterial eq "LI900_A1") { 

$inputf ile«"LI_900_AL_5inputs . in"; 

} 

if <$material eq "LI900_GrEx") ( 

Sinputf ile-"LI_900_GrEx_5inputs . in"; 

) 

if ($material eq "LI2200_A1") { 

$inputf ile-"LI_2200_AL_5inputs. in"; 

) 

if (Smaterial eq "LI2200_GrEx") < 

S input f ile-"LI_2200_GrEx_5inputs. in"; 

} 

if ($material eq "RCC_GrEx") t 

Sinputf ile«"RCC_GrEx_5inputs . in"; 

) 

if ($material eq "RCC_A1"> ( 

Sinputf ile«"RCC_AL_5inputs . in"; 

) 

if (Smaterial eq "SiC_GrEx") { 

Sinputfile-"SIC_GrEx_5inputs .in"; 

> 

if ($material eq "SiC_Al") { 

Sinputf ile«"SIC_AL_5 inputs. in"; 

) 

if ($material eq "TUFI_GrEx"> { 

Sinputf ile-"TUFI_GrEx_5inputs . in"; 

) 

if ($material eq "AFRSI_A1") { 

$input f ile-"AFRSI Al_3inputs . in"; 

} 

if (Smaterial eq "AFRSI_GrEx") { 

Sinputf ile«"AFRSI_GrEx_3inputs . in"; 

> 

if ($material eq "CFBI_A1") { 

Sinputf ile-"CFBI_Al_3inputs . in"; 

) 

if (Smaterial eq "CFBI_GrEx" ) { 
$inputfile«"CFBI_GrEx_3inputs . in"; 

) 

if (Smaterial eq " DURAFRS I_A1 " ) { 

Sinputf ile*"DURAFRSI_Al_3 inputs . in"; 

) 

if (Smaterial eq "DURAFRSI_GrEx" ) { 

Sinputf ile=*" DURAFRS I GrEx_3inputs . in"; 

} 

if (Smaterial eq "PBI_A1"> { 

Sinputf ile*"PBI Al_3inputs . in"; 

) 

if (Smaterial eq "PBI_GrEx") ( 

Sinputf ile-"PBI_GrEx_3inputs . in"; 

} 

if (Smaterial eq "TABIJU") { 

Sinputf ile«"TABI_Al_3inputs . in"; 

> 

if (Smaterial eq "TABI_GrEx" ) { 

Sinputf ile*"TABI_GrEx_3inputs . in"; 

} 

) 

♦print "input file * Sinputf ile\n" ; 

♦ - 

f Execute TCAT and send ADS output into file named "junkoutput" 
# 

# system ("touch inputs. in"); 

♦ system ("chmod 777 inputs. in"); 
system("cp Sinputfile inputs. in"); 

if ($material_option eq 1) { 
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if ($method eq 1) ( 

system (" . /ground. exe") ; 

} else { 

system (" . /ground2 .exe") ; 

I 

) else ( 

if ($material eq "RCC_GrEx") { 

sy3tem("./ground2.exe"> ; ♦ run TCAT 

) 

els if (Smaterial eq "RCC_A1") | 

system (". /ground2 .exe") ; # run TCAT 

} 

elsif (Smaterial eq "SiC_GrEx") ( 

system (" ./ground2 .exe") ; * run TCAT 

) 

elsif <$material eq "SiC_Al") { 

system (". /ground2 .exe") ; # run TCAT 

) 

elsif ($material eq "TUFI^GrEx") { 

system {". /ground2 .exe") ; t run TCAT 

} 

else { 

system (". /ground. exe' ) ; # run TCAT 

1 

) 

# system ("touch $name ( $loop_index]_temphist" ) ; 

# system ("chmod 777 $name l$loop_index)_temphist") ; 
system("mv fort. 90 groundhold_temphist" ) ; 


# system ("touch $name l$loop_index)_heatratehist") ; 

♦ system ("chmod *777 $name [$loop_index]_heatratehist" ) ; 
system("mv fort. 91 groundhold_heatratehist" ) ; 

open (TCATOUTPUT, "<fort . 12") ; # obtain the current thickness of TPS file 
$_*<TCATOUTPUT> ; 

/ A \s* (\SM/; 

$thickness«$l; 
close (TCATOUTPUT) ; 

$ thkns-$ thickness; 

printf FINALOUTPUT "Groundhold - TPS thickness - %6.3f in. \n", $thickness*100/2 . 54 ; 


♦system ("rm point*"); 

# 

f End of heating analysis 

I 

$avg_thickness-$thkns; 

♦for ($loop_index*$start; $loop_index <- Send? $loop_index++) ( 

♦ $avg thickness-$avg thickness+Sthkns [ $loop_index] ; 

♦ ) 

♦$avg_thickness*$avg_thickness/$end; 

open (DEN, "<material density"); # obtain the density of TPS material 
$_-<DEN>; 

/ A \s* ( \ S * ) / ; 

$tps_density-$l? 
close (DEN) ; 

if ( ($sharp eq "yes") £& (Sbodyside eq "windward")) ( 
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$sharp_weight-4 . 6e-3*$bodyarea; #SHARP weight in lbm 
if (Smaterial eq "RCC_GrEx") { 

$tile_unit_weight- (Savg thickness+O . 00254 ) *3 . 2808* ($avg_thickness*Stps_density+0 . 00254*15 

77.8347) / (0.00254+$avg thickness) /O . 4536/3 . 2808/3 . 2808/3 . 2808; #obtain density of tps 

) 

elsif (Smaterial eq "RCC_A1") { 

$tile_unit_weight- ($avg_thickness+0 . 002 54)*3. 2808* ($avg_thickness*$tps_density+0 .00254*15 

77. 8347) /(0.00254+$avg_thickness)/0. 4536/3. 2808/3. 280B/3. 2808; #obtain density of tps 

) 

elsif ($material eq "SiC_GrEx") { 

$tile unit_weight- ($avg_thickness+0 . 00254 ) *3 . 2808* (Savg thickness *$tps_density+0 .00254*24 
00. 007/(0.00254+$avg_thickness)/0. 4536/3. 2808/3. 2808/3. 2808; #obtain density of tps 

) 

elsif (Smaterial eq "SiCJU") ( 

$tile unit weight- (Savg_thickness+0 . 00254) *3.2808* < Savg_thickness*Stps_density+0 . 00254*24 
00. 007/(0.00254+$avg_thickness)/0. 4536/3. 2808/3. 2808/3. 2808; #obtain density of tps 

) 

elsif ($material eq "TUFI_GrEx") { 

$tile_unit weight- ($avg thickness +0 . 00254 ) *3.2808* ($avg_thickness*Stps_density+0 . 00254*13 
13.53)/<0.00254+$avg thickness ) /0 . 4536/ 3 . 280B/3 . 2808/ 3 . 2808; ^obtain density of tps 

) 

else { 

$tile_unit_weight-$avg_thickness*3 . 2808 *$tps_density/0. 4536/3 .2608/3 . 2808/3 . 2808; 
♦obtain density of tps 

} 

$cile_area_to_body_area-Sbcdytile; 
printf FINALOUTPUT "\n"; 

printf FINALOUTPUT "Chosen TPS material was $material\n\n"; 
printf FINALOUTPUT "Nose SHARP TPS: %6.2f lbm\n", $sharp_weight ; 
printf FINALOUTPUT "Smaterial unit weight: 48. 4f lbm/ft2 \n ", $tile_unit_weight ; 
printf FINALOUTPUT "Smaterial TPS area to body area ratio: 44. 2f \n 
", $tile_area_to_body_area; 

} 

if ((Ssharp eq "yes") (Sbodyside eq "leeward")) { 

$sharp_weight-4 . 6e-3*$bodyarea; ISHARP weight in lbm 

Sblanket_unit_weight-$avg_thickness*3 . 2 808* $tps_density/0 .4536/3.2808/3.2806/3,2808; ; 

♦obtain density of tps 
$blanket_area_to_body_area«$bodyblanket; 
printf FINALOUTPUT "\n"; 

printf FINALOUTPUT "Chosen TPS material was $material\n\n" ; 
printf FINALOUTPUT "Nose SHARP TPS: 46. 2f lbm\n", $sharp_weight ; 

printf FINALOUTPUT "Smaterial unit weight: 46. 2f lbm/ft2 \n ", Sblanket_unit_weight; 
printf FINALOUTPUT "Smaterial TPS area to body area ratio: 46. 2f \n 
",$blanket area to body_area; 

) 

if ((Ssharp eq "no") (Sbodyside eq "windward")) ( 
if (Smaterial eq "RCC_GrEx") ( 

Stile unit weight- ($avg_thickness + 0 .00254)*3.2808* (Savg_thickness*Stps_density+0 . 00254*15 

77 .8347) / (0.00254 + $avg thickness ) /0 . 4 536/3 . 2808/3 . 2808/3 . 2808 ; Iobtain density of tps 

) 

elsif (Smaterial eq "RCC_A1") ( 

$tile_unit_weight- (Savg_thickness+0 . 002 54)*3. 2808* ($avg_thickness* $tps_densi ty + 0 . 00254*15 

77.8347) / (0. 002 54+ Savg thickness) /0 . 4 536/3 . 2808/3 . 2 808 /3 . 2808 ; Iobtain density of tps 

) 

elsif (Smaterial eq "SiCJSrEx") ( 

$tile_unit_weight- ($avg_thickness+0 . 00254 ) *3. 280B* (Savg_thickness*Stps_density+0 . 00254*24 
00.00) / (0.00254+$avg thickness) /0 . 4536/3 . 2808/3 . 2808/3 . 2808; Iobtain density of tps 

} 

elsif (Smaterial eq "SiC_Al") { 

$tile_unit_weight- (Savg_thickness+0 .00254)*3.2808* ( $avg_thickness* Stps_density+0 .00254*24 
00. 00)/<0.00254+Savg_thickness>/0. 4536/3. 2808/3. 280B/3. 2808; Iobtain density of tps 
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} 

elsif <$material eq "TUFI_GrEx" ) { 

$til«_unit weight- ($avg_thickness+0. 00254) *3.2808* ($avg_thickness*$tps density*0 . 00254 *13 
13.53) / <0 . 00254+$avg_thickness) /0. 4536/3. 2808/3. 2808/3. 28 08; #obtain density of tps 

) 

else { 

$tile_unit_weight-$avg_thickness*3.2808*$tps_density/0 .4536/3.2808/3.2808/3.2808; 
♦obtain density of tps 

) 

$tile_area_to_body_area-$bodytile; 

print f FINALOUTPUT "Chosen TPS material was $material\n\n" ; 

printf FINALOUTPUT "Smaterial unit weight: %6.2f lbm/ft2 \n ", $tile_unit_weight; 
printf FINALOUTPUT "$material TPS area to body area ratio: %6.2f \n 
",$tile area_to_body_area; 

) 

if ({Ssharp eg "no") It ($bodyside eq "leeward")) ( 

$blanket_unit_weight«$avg_thickness* 3. 2806 *$tps_density/0. 4536/3, 2808/3 . 2808/3 . 2808; ; 

♦obtain density of tps 
$blanket_area_to_body_area-$bodyblanket; 

printf FINALOUTPUT "Chosen TPS material was $material\n\n" ; 

printf FINALOUTPUT "Smaterial unit weight: %6.2f lbm/ft2 \n ", Sblanket_unit_weight; 
printf FINALOUTPUT "Smaterial TPS area to body area ratio: %6.2f \n 
", $blanket_area_to_body_area; 

) 

print "average thickness: $avg_thickness\n* ; 
print "TPS density: $tps_density\n\n"; 

if <$material_l eq "RCC_Laninate_Surf ace" ) { 


$tile_unit_weight- ($avg_thickness+0 .00254) *3.2808* ($avg_thickness*$tps_density+0 . 00254*15 
77.8347)/ (0. 00254 + $avg_thickness ) /0. 4536/3. 2808/3. 2808/3. 2808; ^obtain density of tps 

} elsif ($material_l eq "SiC_Laminate_Sur f ace" ) i 

$tile_unit_weight- ($avg_thickness+0 . 00254 ) *3.2808* ($avg_thickness*$tps_density+0 .00254*24 
00. 00)/{0.00254+$avg_thickness)/0. 4536/3. 2808/3. 2B08/3. 2808; ♦obtain density of tps 

) elsif ($material_l eq "TUFI_Composite_Layer") { 

$tile_unit_weight« {$avg_thickness+0 . 00254) *3.2808* ($avg_thickness*$tps_density+0 . 00254 *13 
13. 53)/(0.00254+$avg_thickness)/0. 4536/3. 2808/3. 2808/3. 2808; ^obtain density of tps 

) else ( 

$t ile_unit_weight-$ avg_thickne ss* 3 . 2808 *$tps_density/0 . 4536/3 . 2808/3 .2808/3.2808; 
♦obtain density of tps 

) 

printf FINALOUTPUT "\n"; 
if ($method eq 1) { 

printf FINALOUTPUT "Chosen TPS material was $material_l\n\n"; 

} else { 

printf FINALOUTPUT "Chosen TPS material was $material_2\n\n"; 

) 

printf FINALOUTPUT "TPS unit weight: %8.4f lbm/ft2 \n ", $tile_unit_weight; 
close (FINALOUTPUT) ; 
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♦ 

# print needed information to screen and file named "thickness.txt" 

I 

open (OUT, "<thickness. txt") ; #opens the input file for TCAT script 

open (WEB, "<tcat_output . html") ; #open teat html file 

♦print "Content-type : text/html\n\n"; 
while (<WEB>) { 

if {/Insert stuff here/) { 
print "<BR>\n"; 
while (<OUT>) ( 

print "<BR>\n"; 
print $_; 

} 

) 

else{ 

print $_; 

) 

} 

close (WEB); 
close (OUT) ; 

# ) 
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Appendix I: ‘compare.cgT Ground Hold Multiple Analysis CGI Script 

# ! /usr/local/bin/perl 

I - 

# Obtain Information From Web Browse 

# - 

# This section of code allows a cgi script to run from a web browser. 

# The information is passed using the POST and the input information is 

# parceled out in a mesh type variable and stores in $FORM( $name } . 

I 

$PATH« 'pwd' ; 
chop {$ PATH) ; 

read (STDIN, $buffer, $ENV{ ’ CONTENT_LENGTH * ) ) ; 

Gpairs - split ( / * / r $buffer ) ; 
foreach $pair (Spairs) 

I 

($name r $value) - split (/-/, $pair) ; 

$value — tr/ +/ /; 

$value s/% Ua-fA-FO-9) (a-fA-FO-9) ) /pack ("C", hex($l>)/eg; 

$value — s /-!/ -!/g; 

$FORM{$name) - $value; 

) 

# - — 

# assign information from input website to appropriate variable 

# 

$material_option*$FORM{rr.aterial_option) ; 

$matnum[l]-$FORM(matnum_l) ; 

$matnuml2]-$FORM{matnum_2) ; 

$matnum(3 J«$FORM(matnum_3) ; 

$material_l [l)-$FORM(material_l_l ) ; 

$material_2 [ 1 ) -$ FORM (mater ial_2_l ) ; 

$material_3 [ 1 ] «$FORM (material_3_l ) ; 

$material_4 [ 1 ] “$FORM (material_4_l ) ; 

$material_5 [ 1 ] -SFORM (material_5_l ) ; 

$material_l [2)«$F0RM{material_l_2 ) ; 

$material_2 [2)*$FORM{material_2_2 ) ; 

$material_3 [ 2 ] -$FORM{ materia 1_3_2 } ; 

$material_4 [2 ] =$FORM{ materia 1_4_2 J ; 

$material_5 [2 ] =$ FORM { mater ial_5_2 ) ; 

$material_l (3]=$FORM(material_l_3); 

$material_2 [3)-$FORM (mater ial_2_3 } ; 

$material_3 ( 3] -$FORM{ materia 1_3_3 } ; 

$material_4 [ 3)-$FORM(material_4_3) ; 

$mater ial_5 I 3 ) -$FORM (material_5_3 ) ; 

Smethod I 1 ) *$FORM (method_l } ? 

$method [ 2 ] -$FORM { method_2 ) ; 

$met hod ( 3 ] - $ FORM ( me t hod_3 ) ; 

$holdtime=$FORM ( holdtime } ; 

$t ime step- $ FORM ( time step ) ; 

$ambient-$FORM{ ambient ); 

$fueltemp=$FORM{ fuel temp) ; 

$frost«$FORM{ frost ) ; 

open (TIME, ">time") ; 

print TIME *$holdtime ! /Time of Groundhold"; 

close (TIME ) ; 

Sheatrate - 1; 
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Sheattransfer - $heatrate • $ambient *3.15459074506; 
open (HEATXFER, ">heat ") ; 

print HEATXFER "Sheattransfer !/Heat Flux"; 

close (HEATXFER) ; 

$ambient - ($ambient-32 ) * (5/9) +273 . 15; 

$f rost - ($frost-32)*(5/9)+273.15; 

$f ueltemp - (Sfueltemp-32 ) • (5/9) +273 . 15; 
open (TANK, ">fueltemp") ; 

print TANK "Sfueltemp ’/Temperature of fuel"; 

close (TANK) ; 


for ($case«l; $case <- 3? $case++) ( 
if ($material_option eq 1) ( 

for ($loopnum-l; $loopnum <« $matnum { $case ] ; $loopnum++) ( 
if ($loopnum eq 1) ( 

open (DAT, "<$material_l [$case] ") ; 

$i - 0; 

while (<DAT>) { 

Si - Si + 1; 
my(Sline) ■ $_; 
chomp (Sline) ; 

if ($i eq 1) ( 

Snodestl] » Sline; 

} elsif (Si eq 2) ( 

Scond ( 1 ] - Sline; 

} elsif (Si eq 3) { 

Sdensity(l) - Sline; 

) elsif (Si eq 4) { 

$specheat[l] - Sline; 

) elsif (Si eq 5) { 

Sthick [ 1 ] * Sline; 

) elsif (Si eq 6) { 

$ temp limit ( 1 ] -Sline; 

) elsif (Si eq 7) ( 

Semis sivity- Sline; 

) 

} 

close (DAT) ; 

) elsif (Sloopnum eq 2) { 

open (DAT, "<$mater ial_2 [Scase] ") ; 

Si - 0; 

while (<DAT>) ( 

Si «■ Si + 1; 
my (Sline) - $_; 
chomp ($line) ; 

if (Si eq 1) ( 

$nodes[2] - Sline; 

} elsif (Si eq 2) { 

$cond[2] - Sline; 

} elsif (Si eq 3) { 

$density{2] - Sline; 

} elsif (Si eq 4) ( 
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$specheat[2] ■ $line; 
) elsif ($i eq 5) ( 
$thick[2] - $line; 

} elsif ($i eq 6) { 
$templimit [2]«$line; 

} 


close (DAT) ; 

) elsif ($loopnum eq 3) { 

open (DAT, "<$material_3 IScase] ") ; 
$i - 0 ; 

while (<DAT>) { 

$i - $i + 1; 
my($line) - $_; 
chomp ($line) ; 

if ($i eq 1) ( 

$nodes(3] - $line; 

) elsif ($i eq 2) { 

$cond[3] - $line; 

) elsif ($i eq 3) i 
$density(3] * Sline; 

) elsif ($i eq 4) { 

$specheat(3J - $line; 

} elsif ($i eq 5) ( 

$thick [ 3 ) - $line; 

) elsif ($i eq 6) { 
Stemplimit [3]«$line; 


close (DAT) ; 

} elsif ($loopnum eq 4) ( 

open (DAT, w <$material_4 [$casej ") ; 
$i - 0; 

while (<DAT>) { 

$i - $i + 1; 
my(51ine) - 
chomp ($line) ; 

if ($i eq 1) { 

$nodes[4J - $line; 

) elsif <$i eq 2) { 

$cond[4] - $line; 

) elsif ($i eq 3) { 
$density[4] * $line; 

) elsif ($i eq 4) ( 
$specheat[4J - $line; 

) elsif ($i eq 5) { 

$thick [4 ] - $line; 

} elsif ($i eq 6) { 

$ temp limit [4 ) -$line; 


close (DAT) ; 

} elsif (Sloopnum eq 5) { 
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open (DAT, "<$material_5 [ Scase] ") ; 
$i - 0; 

while (<DAT>) ( 

Si « Si + 1; 
my ($line) - $_; 
chomp ($line) ; 

if ($i eq 1) { 

$nodes[5] - $line; 

) elsif (Si eq 2) ( 

$cond[5) * Sline; 

) elsif (Si eq 3) ( 
$density[5] - Sline; 

} elsif (Si eq 4) ( 
$specheat(5] - Sline; 

) elsif (Si eq 5) ( 

$thick{ 5] - Sline; 

) elsif (Si eq 6) { 
Stemplimit [5] -Sline; 


close (DAT) ; 


) 


) 

#########♦ 

# Create the TPS stackup file by reading the properties of each material selected and 
entering # 

# their values into a input file recognizable by T-CAT. 

# 

*###*###«* 


# system("touch stackup . in" ) ; 

# system ("chmod 777 stackup. in" ) ; 
open (STACKUP, ">stackup. in") ; 
print STACKUP "Smatnum ( Scase ] 
print STACKUP "Stimestep 


♦create TPS stackup file 

♦change permissions of TPS stackup file 

! /number of materialsXn" ; 

I/time step\n"; 


for ($loopnum*l; Sloopnum <« Smatnum [Scase) ; $loopnum4+) ( 

printf STACKUP "%i ! /number of nodes of material 

$loopnum\n", Snodes [Sloopnum) ; 

> 

print STACKUP "Semissivity !/emissivity of material l\n"; 

for <$loopnum*l; Sloopnum <- Smatnum [Scase) ; $loopnum++) { 

printf STACKUP "%f [/thermal conductivity of material 

$loopnum\n", Scond [Sloopnum] ; 

) 

for ($loopnum«l; Sloopnum <- Smatnum [ Scase) ; $loopnum44) ( 

printf STACKUP "%f {/density of material $loopnum\n", Sdensity [ Sloopnum] ; 

) 

for (Sloopnum-1; Sloopnum <- $matnum[Scase) ; Sloopnum+4) ( 
printf STACKUP "%f ! /specific heat of material 

SloopnumVn", Sspecheat [Sloopnum] ; 

) 

print STACKUP "Sambient [/initial temp (deg K)\n"; 

for (Sloopnum-1; Sloopnum <* Smatnum [Scase] ; Sloopnum44) { 


Stephen James Izon 


Pace 52 of 81 



Improvements in Thermal Protection 


L»£ Ca pab ilities 


Systems jfoi 


printf STACKUP "%f ’/TPS thickness of material $loopnum\n" , $thick [ $loopnum] ; 

) 

print STACKUP "$frost ! /Temperature limit of material lXn"; 

print STACKUP "Stemplimit [ 1 ] I /Temperature limit of material l\n"; 

for ($loopnum«2; $loopnum <- $matnum l$case] -1; $loopnum++) ( 
printf STACKUP "%f ’/Temperature limit of material 

$loopnumXn", $templimit [$loopnum] ; 

printf STACKUP "%f ! /Temperature limit of material 

$loopnum\n", $templimit l$loopnum] ; 

) 


print STACKUP 
SmatnumXn"; 

print STACKUP 
print STACKUP 
print STACKUP 
print STACKUP 
print STACKUP 


"$templimit [$matnum($case] ] ’/Temperature limit of material 

"$fueltemp !/backface temperature liraitXn"; 

"6 !/istrat parameterXn"; 

*5 !/iopt parameterXn"; 

"5 ! /ioned parameterXn" ; 

"1340 !/iprint parameterXn"; 


close (STACKUP) ; 

) 


############ 

# This section of the program creates a mock MINIVER file. The three columns of data 
created are I 

I the time, ambient temperature, and rate of heat transfer. 

I 

########################################M############################################### 
III Ml###### 

$numtimes « $holdtime/$timestep + 1; 

$time I 1 ] - 0; 

for ($i“2; $i <- Snumtimes; $i++) { 

$time[$i] « $time[$i-l] + $timestep; 


) 

#system("touch l_stackup.s") ; #create mock MINIVER file 

#system("chmod 111 l_stackup.s") ; Ichange permissions of file 

open (MINI, ">l_stackup . s") ; 

print MINI " Mock MINIVER file for input into ground hold analysis 10 

l.OXn"; 

for ($i-l; $i <» $numtimes; $i++) ( 

print MINI "$time($i] $ambient $heattransfer\n"; 


) 

print MINI "-100. Un- 
close (MINI); 

# — 

# MINIVER file is partitioned into body point files 

# partition the MINIVER file into separate body point files that include time, 

# heat rate, and radiation equilibrium temperature 

# 

$i-0; Initialization of counter index "i" 

$f ilecount-0; 

open (FILE, "<l_stackup. s") ; #open statement with filehandle FILE 

until (eof FILE) { Igoes through the MINVER file until the end is 

reached 
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START: $_-<FILE>; fsk ips the first line of the data file 

if (/* . <73 )(VS*> /) ( fsearches for the bodypoint label 

$title-$l; 

# system ("touch groundhold" ) ; fcreate body point file 

# system ("chmod 777 groundhold"); #change permissions inorder to write to body point 

file 

$f ilecount-$f ilecount+1; 

) # end of if statement 

while (<FILE>) { #while loop that reads thru each line of the file 

$i~$i+l ; 

/'Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS*) Vs* (VS* 
) Vs* (VS* ) Vs* (VS*)/; 

$time ( $i] -$1; #array declaration for time 

Stradeq ($i) -$2; f array declaration for rad eq temp 

$qconv($i]-$3; farray declaration for conv heat rate 

$n-$i; 

if ($1 — 100) ( 

open (NEWFILE, ">groundhold" ) ; #open statement for new file "newsample" 
$timesteps-$n-l ; 

printf NEWFILE "$timesteps Vn"; #formatted print statement" 

f or < $i*l ; $i<$n; $i++) { 

printf NEWFILE "If Vn" , $time [ $i] ; 
printf NEWFILE "If Vn" , $qconv { $i } ; 
printf NEWFILE "If Vn" , $t radeq [ $i ] ; 

) #end of for loop 

close (NEWFILE) ; 

$i«0; 

goto START 

) lend of if statement 

} lend of while loop 
) fend of until loop 
close (FILE); 

# — 

# incrementally move each body point file into "miniver. in" and execute teat 

t perl script that creates the filenames needed inorder to batch execute teat 
t filename are created via input from the user 

# inputs (file start_number , f ile_end_number, body side) 

t — : : : 

$start*l; tstart of the file index 

#$end« ($f ilecount-1) /2 ; fend of the file index 

$end*l ; 

open (FILE, ">inputs_f or_output s" ) ; 
printf FILE "3Vn"; 
close (FILE) ; 

fsystem ( "touch thicfcness.txt"); 

♦system ( "chmod 777 thickness.txt"); 
open (FINALOUTPUT, ">thickness . txt" ) ; 

print FINALOUTPUT "TPS Ground hold design resultsVnVn"; 

# 

# system ( "touch miniver . in") ; 

f system ( "chmod 777 miniver . in") ; 
system("cp groundhold miniver . in") ; 

# 

# assign material properties 

# 

if ($material_option eq 1) ( 
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Sinputf ile - "stackup. in**; 

) else { 

if {$material eq "AETB8_Al"> { 

Sinputf ile-"A£TB_8_AL_5inputs . in"; 

) 

if (Smaterial eq "AETB8_GrEx") { 

$inputf ile-"AETB_8_GrEx_5 inputs . in"; 

) 

if (Smaterial eq "AETB12_A1" ) { 

Sinputf ile-"AETB_12 AL_5inputs . in"; 

) 

if (Smaterial eq "AETB12_GrEx" ) ( 

Sinputf ile«"AETB_12 GrEX_5inputs . in"; 

} 

if ($material eq "FRCI12_A1") { 

Sinputf ile-"FRCI_12_AL_5inputs . i"; 

) 

if (Smaterial eq "FRCI12_GrEx") ( 

$inputf ile-"FRCI_12_GrEx_5inputs . in" ; 

) 

if (Smaterial eq "FRCI20JU") { 

Sinputf ile-" FRCI_20_AL_5 input s . in" ; 

) 

if (Smaterial eq "FRCI20_GrEx" ) { 

$ input file*" FRCI_20_GrEx_5 inputs - in* ; 

) 

if ($material eq "LI900_A1") { 

$inputf ile-"LI_900_AL 5inputs . in" ; 

) 

if (Smaterial eq "LI 900_GrEx") { 

Sinputf ile- "LI_900_GrEx_5 inputs . in" ; 

} 

if (Smaterial eq "LI2200_A1") < 

Sinputf ile- "LI_2 2 00_AL 5 inputs .in"; 

} 

if (Smaterial eq "LI2200_GrEx" ) { 

Sinputf i le- " LI_2 2 00_GrEx_5 inputs .in"; 

} 

if (Smaterial eq "RCC_GrEx") { 

Sinputf ile="RCC_GrEx_5 inputs . in" ; 

1 

if (Smaterial eq "RCC_A1") ( 

Sinputf ile="RCC_AL_5 inputs .in"; 

} 

if (Smaterial eq "SiC_GrEx") { 

Sinputf ile«"SIC_GrEx_5inputs . in"; 

) 

if (Smaterial eq "SiC_Al") { 

Sinput file-"SIC_AL_5inputs .in"; 

) 

if (Smaterial eq "TUFI_GrEx") { 

Sinputf ile- "TUFI_GrEx_5 inputs . in"; 

) 

if (Smaterial eq "AFRSI_A1") ( 

Sinputf ile-"AFRSI_Al_3inputs . in"; 

} 

if (Smaterial eq "AFRSI_GrEx" ) ( 

Sinputf ile«"AFRSI_GrEx_3inputs , in"; 

} 

if (Smaterial eq "CFBI_A1"> { 

Sinputf ile- "CFBI_A1_3 inputs . in"; 

} 

if (Smaterial eq "CFBI_GrEx") { 
Sinputf ile«"CFBI_GrEx 3inputs . in"; 

) 

if (Smaterial eq "DURAFRSI_A1") ( 
Sinput f i le-"DURAFRSI_Al_3 inputs . in"; 

) 

if (Smaterial eq "DURAFRSI_GrEx" ) ( 
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$inputf ile«"DURAFRSI__GrEx_3 inputs .in"; 

) 

if ($mater ial eq "PBI_A1") { 

$inputf ile«"PBI_Al_3 inputs . in" ; 

) 

if ($material eq "PBI_GrEx"> { 

$inputfile-"PBI_GrEx 3inputs . in"; 

) 

if ($material eq "TABIJU") { 

$inputf ile«" TAB I_A1_3 inputs . in" ; 

) 

if (Smaterial eq "TABI_GrEx") { 

$inputf ile«"TABI GrEx_3inputs . in" ; 

} 

) 

* - 

# Execute TCAT and send ADS output into file named "junkoutput" 

# - 


I system("touch inputs. in"); 

# system ("chmod 777 inputs. in"); 

system("cp $inputfile inputs. in » error 2>4l")? 


# system ("touch inputs [ $case] "} ; 

# system ("chmod 777 inputs ( $case] ") ; 
system ("cp inputs. in inputs ($casej ") ; 

if ($material_option eq 1) ( 

if ($method($case) eq 1) ( 

system (" . /ground. exe") ; 

) else { 

system (" . /ground2 .exe") ; 

) 

) else ( 


> 


if ($material eq "RCC_GrEx") ( 
system("go4 ,exe>junkoutput") ; 

) 

elsif ($material eq "RCC_A1") ( 
system ("go4 .exe>junkoutput" ) ; 

) 

elsif ($material eq "SiC_GrEx") ( 
system ("go4 ,exe> junkoutput" ) ; 

) 

elsif ( $material eq "SiC_Al") { 
system ( "go4 . exe> junkoutput" ) ; 

) 

elsif ($material eq "TUFI_GrEx") ( 
system ("go4 . exe> junkoutput") ; 

) 

else ( 

system ( "go3 . exe> junkoutput") ; 

) 


# run TCAT 

# run TCAT 
I run TCAT 

# run TCAT 

# run TCAT 
I run TCAT 


# system ("touch $name l$loop_index]_temphist") ; 

I system ("chmod 777 $name ( $loop_index]_temphist" ) ; 

system ("mv fort. 90 groundhold_temphist ( $case] >> error3"); 


# system("touch $name [$loop_index]_heatratehist ") } 

# system ("chmod 777 $name ($loop__index]_heatratehist") ; 
system ("mv fort. 91 groundhold_heatratehist [Scase] >> error4"); 
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open (TCATOUTPUT, "<fort . 12") ; # obtain the current thickness of TPS file 
$_«<TCATOUTPUT>; 

/ A \s* ( \S* ) /; 

$thickness*$l ; 
close (TCATOUTPUT) ; 

$thkns-$thickness; 

printf FINALOUTPUT "Groundhold - TPS thickness - %6.3f in . \n" , $thickness*100/2 . 54; 


#s ystem("rm point*"); 

4 — 

4 End of heating analysis 

4 - 

$avg_thickness-$thkns; 

4for ($loop_index-$start; $loop_index <- $end? $loop_index++) ( 

# $avg thickness«$avg thickness+$thkns ( $loop_indexT; 

#} 

#$avg_thickness-$avg_thickness/$end; 

open (DEN, "<material_density") ; # obtain the density of TPS material 
$_«<DEN>; 

/ A \s* (\S* ) /; 

$tps_density-$l; 
close (DEN) ; 

if (($sharp eq "yes") ti (Sbodyside eq "windward")) ( 

$sharp_veight-4 . 6e- 3*$bodyarea; ISHARP weight in lbm 
if ( $ma terial eq "RCC_GrEx") ( 

$tile_unit_weight- ($avg_thickness+0 . 00254 > *3 .280B* ($avg_thickness*$tps_density+0 .00254*15 
77.8347) / (0. 00254+ $avg_thickness) /0. 4536/3. 2808/3. 2808/3. 2808; Iobtain density of tps 

) 

elsif ($material eq "RCC_A1") { 

$tile_unit_weight- ($avg_thickness+0 . 00254) *3.2808* ($avg_thickness*$tps_density+0 .00254*15 
77.8347)/ (0. 00254+$avg_thickness) /0. 4536/3. 2808/3. 2808/3. 28 08; Iobtain density of tps 

} 

elsif ($material eq "SiC_GrEx") ( 

$tile_unit_weight- ($avg_thickness+0 . 00254) *3.2808* ($avg_thickness*$tps_density+0 . 00254 *24 
00 . 00) / (0 . 00254+$avg_thickness) /0 . 4536/3 . 2808/3 . 2808/3 . 2808; Iobtain density of tps 

1 

elsif ($material eq "SiC_Al") ( 

$tile_unit_weight- ($avg_thickness40 . 00254) *3.2808* ($avg_thickness*$tps_density+0 .00254*24 
00. 00)/(0.00254+$avg_thickness)/0. 4536/3. 2808/3. 2808/3. 2808; Iobtain density of tps 

) 

elsif ($material eq **TUFI_GrEx" ) ( 

$tile_unit_weight- ($avg_thickness+0 . 00254) *3.2808* ($avg_thickness*$tps_density+0 .00254*13 
13. 53)/(0.00254+$avg_thickness>/0. 4536/3. 2808/3. 2808/3. 2B08; Iobtain density of tps 

> 

else { 

$tile_unit_weight-$avg_thickness*3 . 2808*$tps_density/0 . 4536/3 . 2808/3 .2808/3.2808; 
♦obtain density of tps 

) 

$tile_area_to_body_area*$bodytile; 
printf FINALOUTPUT "\n"; 

printf FINALOUTPUT "Chosen TPS material was $material\n\n"; 
printf FINALOUTPUT "Nose SHARP TPS: %6.2f lbm\n", $sharp_weight; 
printf FINALOUTPUT "$material unit weight: %8.4f lbm/ft2 \n ", $tile_unit_weight; 
printf FINALOUTPUT "$material TPS area to body area ratio: %4.2f \n 
", $tile_area_to_body_area; 


) 

if ( (Ssharp eq "yes") && ($bodyside eq "leeward")) { 

$sharp_weight«4 . 6e-3*$bodyarea; ISHARP weight in lbm 
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$blanket_unit_weight-$avg_thickness*3. 2808*$tps_density/0 . 4536/3 . 2808/3 .2808/3 .2808 ; ; 

♦obtain density of tps 
$blanket_area_to_body_area«$bodyblanket; 
printf FINALOUTPUT "\n"; 

print f FINALOUTPUT "Chosen TPS material was $material\n\n" ; 
printf FINALOUTPUT "Nose SHARP TPS: %6.2f lbm\n", $sharp_weight ; 

printf FINALOUTPUT "$material unit weight: %6.2f lbm/ft2 \n ", $blanket_unit_weight ; 
printf FINALOUTPUT "$material TPS area to body area ratio: %6.2f \n 
", $blanket_area to_body_area; 

) 

if { ($sharp eq "no") && ($bodyside eq "windward")) 1 
if <$material eq "RCC_GrEx") < 

$tile_unit_weight- ($avg_thickness+0 .00254)*3.2808* ($avg_thickness*$tps_density+0 .00254*15 

77.8347) / (0 . 00254 + $avg^thickness) /0. 4536/3. 2808/3. 2808/3. 2808; ♦obtain density of tps 

> 

elsif ($material eq "RCC_A1") { 

$tile_unit_weight« <$avg thickness+0 . 00254) *3 . 2808* {$avg_thickness*$tps_density+0 .00254*15 

77.8347) / (0.002S4 + $avg thickness) /0 . 4 536/3 . 2808/3 . 2806/3 . 2808 ; lobtain density of tps 

> 

elsif ($material eq "SiC_GrEx") ( 

$tile_unit_weight« ($avg_thickness+0 .00254)*3.2808* ($avg_thickness*$tps_density+0 .00254*24 
00. 00)/{0.00254+$avg_thickness)/0. 4536/3. 2808/3. 2808/3. 2808; #obtain density of tps 

) 

elsif (Smaterial eq "SiC_Al") { 

Stile unit weight^ ($avg_thickness+0 .00254)*3.2808* ( $avg_t hickness * $tps_dens it y+0 .00254*24 
CO. 00) / (0.00254+$avg_thickness) /0 . 4536/3. 2808/3. 2808/3. 2808; #obtain density of tps 

} 

elsif (Smaterial eq "TUFi_GrEx") { 

Stile unit_weight* ($avg_thickness-*-0 . 002 54 )*3. 2808* ($avg_thickness*$tps_density+0 . 002 54*13 
13.537/ (0. 00254+$avg_thickness) /0. 4536/3. 2808/3. 2808/3. 2808; #obtain density of tps 

} 

else { 

$tile_unit_weight*$avg_t hickness* 3 .2808*$tps_density/0 .4536/3.2808/3.2806/3.2808; 
♦obtain density of tps 

) 

$tile_area_to_body_area~$bodytile; 

printf FINALOUTPUT "Chosen TPS material was $mater ial\n\n" ? 

printf FINALOUTPUT "$material unit weight: %6.2f lbm/ft2 \n " , $tile_unit_weight ; 
printf FINALOUTPUT "$material TPS area to body area ratio: %6.2f \n 
" , $tile_area_to_body_area; 

) 

if ( ($sharp eq "no") ll ($bodyside eq "leeward")) { 

$blanket_unit_weight«$avg_thickness*3 . 2808*$tps_density/0 .4536/3.2808/3.2808/3.2808;; 

♦obtain density of tps 
$blanket_area_to_body_area-$bodyblanket ; 

printf FINALOUTPUT "Chosen TPS material was $material\n\n" ; 

printf FINALOUTPUT "$material unit weight: %6.2f lbm/ft2 \n ", $blanket_unit_weight ; 
printf FINALOUTPUT "$material TPS area to body area ratio: %6.2f \n 
" , $blanket_area_to_body area; 

) 

♦print "average thickness: $avg_thickness\n"; 

♦print "TPS density: $tps_density\n\n"; 


if ($material_l [$case] eq "RCC_Laminate_Surface" ) { 

$tile_unit_weight« ($avg_thickness+0 . 00254) *3.2808* ($avg_thickness* $tps_density+0 .00254*15 

77.8347)/ (0. 00254 + $avg_thickness) /0. 4536/3. 2808/3. 2808/3. 2808; lobtain density of tps 

) elsif ($material_l [$casej eq "SiC_Laminate_Surface") ( 


$tile_unit_weight« ($avg_thickness+0 . 002 54 )*3. 2808* ($avg_thickness*$tps_density+0 .00254*24 
00. 00)/(0.00254+$avg_thickness)/0. 4536/3. 2808/3. 2B08/3. 2808; lobtain density of tps 
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) elsif ($material_l (Scase] eq "TUFI_Composite_Layer") { 


$tile_unit_weight« (Savg thickness+O . 00254) *3.2808* ($avg_thickness*$tps density+0 . 00254 *13 
13. 53>/(0.00254+$avg_thickness}/0. 4536/3. 2B08/3. 2808/3. 2808; #obtain density of tps 

) else ( 

$tile_unit_weight*$avg_thickness*3 . 2808*$tps_density/0 .4536/3.2808/3.2808/3.2808; 
•obtain density of tps 

} 

printf FINALOUTPUT "\n"; 
if ($method[$case] eq 1) { 

printf FINALOUTPUT "Chosen TPS material was $material_l [ Scase) \n\n"; 

) else { 

printf FINALOUTPUT "Chosen TPS material was $material_2 (Scase) \n\n"; 


) 

printf FINALOUTPUT "TPS unit weight: %8.4f lbm/ft2 \n ", $tile_unit_weight ; 


close (FINALOUTPJT) ; 

$ race thick [ $case] "$avg_thickness; 

$casedensity [ Scase] *$tps_density; 

$unitweight [ $case] «$tile_unit_weight; 

system{"mv thickness . txt thickness [$case] . txt >> error4") ; 


) 

if (< Sunitweight ( 1 ] 
$min « 3; 

) 

if ( (Sunitweight 1 1 } 
$min ■ 2; 

} 

if ( ( Sunitweight [ 2 ] 
$min • 3; 

} 

if (( Sunitweight [2] 
Smin - 1; 

if (( Sunitweight [ 3} 
Smin - 1; 

) 

if ( (Sunitweight [ 3] 
Smin « 2; 

) 


gt Sunitweight [2] ) 
gt Sunitweight ( 3 J ) 
gt Sunitweight ( 1 ] ) 
gt Sunitweight (3) ) 
gt Sunitweight [ 2 ] ) 
gt Sunitweight (1 ] ) 


&& (Sunitweight [2 J 
&& ( Sunitweight [ 3] 
(& ( Sunitweight [ 1 ] 
< Sunitweight [ 3] 
&& (Sunitweight [2] 


gt Sunitweight 13] ) ) ( 
gt Sunitweight [2 ]) ) { 
gt Sunitweight [ 3] ) ) { 
gt Sunitweight [1 ]) ) { 
gt Sunitweight [1 ]) ) { 


&& (Sunitweight ( 1] gt Sunitweight [2 ]) ) ( 


open (RESULTS, ">results . txt" ) ; 

printf RESULTS "The optimum design for this analysis is configuration $min\n"; 
printf RESULTS "The configuration consists of the following materials : \n\n"; 

printf RESULTS " Material 1: $material_l ( Smin] \n"; 

if (Smatnum [ Smin] gt 1) { 

printf RESULTS " Material 2: $material_2 l$min] \n"; 

} 

if (Smatnuml Smin] gt 2) { 

printf RESULTS " Material 3: $material_3 [Smin] \n"; 
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) 

if ($matnuml $min] gt 3) { 

printf RESULTS • Material 4: $material_4 [$min] \n"; 

} 

if (Smatnum [ $minj gt 4) { 

printf RESULTS " Material 5: $material_5 [$min] \n"; 

) 

$finalweight«$casethick!$min] *3. 2808*$casedensity($min]/0. 4536/3. 2B08/3. 2808/3. 2808; 
printf RESULTS "\n"; 

printf RESULTS " TPS Ground hold design results \n\n"; 

printf RESULTS "You chose to optimize layer $method [$min] of the TPS stack. \n\n"; 
if ($method[ Smin] eq 1) { 

printf RESULTS "The TPS material for layer Smethod [$min] is $material_l { Smin] . \n\n" ; 
printf RESULTS "$material_l l $min] thickness - 16. 3f 
in . \n\n" f $case thick! Smin] *100/2.54; 

printf RESULTS "TPS unit weight: %8.4f lbm/ft2 \n ", Sfinalweight ; 

) else { 

printf RESULTS "The TPS material for layer Smethod [$min] is $material_2 ( $min] . \n\n" ; 
printf RESULTS "$mater ial_2 [ Smin] thickness - %6.3f 
in . \n\n" , $casethick [ $min] *100/2.54; 

printf RESULTS "TPS unit weight: 45. 4f lbm/ft2 \n " , $f inalweigh ; ; 

} 

close (RESULTS) ; 

# - 

# print needed information to screen and file named "thickness.txt" 

I - 

open (OUT, "<results . txt" ) ; lopens the input file for TCAT script 

open (WEB, "<t cat_output . html" ) ; lopen teat html file 

print "Con tent- 1 ype : text /html \n\n" ; 
while (<WEB>) ( 

if (/Insert stuff here/) ( 
print "<BR>\n"; 
while (<OUT>) ( 

print "<BR>\n"; 
print 

) 

) 

else ( 
print 

} 

} 

close (WEB ) ; 
close (OUT) ; 

» ) 
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Appendix J: ‘ground.F FORTRAN Source Code 

program ground 


• Program to solve the one dimensional heat equation through a stackup 

• of disparate TPS materials. This program finds the optimum thickness • 

• of the surface layer of the TPS material to reduce the TPS unit • 

• weight, while eliminating frosting at the surface of the TPS stackup. • 


implicit none 

integer i, j, n, nodes (6) , suml,sum2, sum3, totnodes,matnum, count 
integer iter,maxiter 

real *8 dt, k (6) ,den (6) , cp<6) , initT,L(6) , dx (6) , sigma, epsi,eps 
1, alpha {6) , lam (6) , a (400) , b (400) ,c (400) ,T (10000, 400) ,noraf 
l,normdelT, f (400) ,delT(400) , twonorm, time (10000) 
l,qcond (400) , qrad(400) ,qcv(400) , g (24) , obj , x (12) ,dTemp 
1, kair,Tb(1000) , error, diff (1000) ,templmt (10) ,backT,heat 
1, root (1000) , toler 


• Obtains information from the file 'inputs, in’ which contains TPS • 

• material properties, used in the solution of the heat equation. • 


open {ID, file-' inputs. in*) 
read (10,*) matnum 
read(10,*) dt 
do i-1, matnum 
read(10,*) nodes(i) 
enddo 

read(10,*) epsi 
do i-1, matnum 
readdO,*} k(i) 
enddo 

do i-1, matnum 
read (10, •) den(i) 
enddo 

do i-1, matnum 
read{10,*) cp(i) 
enddo 

readdO,*) initT 
do i-1, matnum 
readdO,*) L(i) 
enddo 

do i»l,2*matnum 
read(10,*) templmt(i) 
enddo 


(/number of materials 
!/time step 

(/number of nodes per material 
! /emmissivity of surface 
(/thermal conductivities 

! /densities 

! /specific heats 
(/initial temperatures 
! /thicknesses 

(/temp limits for each material 


• Writes the density of the surface TPS material for use in calculating 

• the TPS material unit weight in the perl script. • 

• Also writes the heat transfer rate and the backface fuel temperature • 

• to files for use in the one-dimensional heat equation solver. 


open (40, file-’ heat ' ) 
read (40, * ) heat 
close (40) 

open (41, f ile»* fueltemp' ) 
read (41,*) backT 
close (41) 

open (9, file«*material_density* ) 
write ( 9, *) den(l) 
close (9) 
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* Initializes optimization parameters. 


toler-0.0001 
error-1 . 
maxiter-100 


* Solves the heat equation for using two initial guesses for the * 

* thickness of the TPS. Values obtained for the objective function * 

* will be sent to the next while loop which implements the * 

* Newton-Raphson method to converge to the solution. * 


if (matnum.eq. 1) then 

call one (matnum, dt, nodes, epsi, k, den, cp, initT,L,T,n, 

+ time,qcv, qcond, qrad, backT, heat) 

el seif (matnum.eq. 2) then 

call two (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum .eq. 3) then 

call three (matnum, dt, nodes, epsi, k, den, cp, initT, L,T, n, 
+ time, qcv, qcond, qrad, backT, heat ) 

elseif (matnum.eq. 4) then 

call four (matnum, dt., nodes, epsi, k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat ) 

elseif (matnum.eq. 5) then 

call five (matnum, dt, nodes, epsi, k, den, cp, initT, L,T, n, 

♦ time, qcv, qcond, qrad, backT, heat) 

endif 

diff (1) -T (n, 1 ) -templmt (1) 
root (1 ) -L (1) 

L<1)«0. 00001 

if (matnum.eq. 1) then 

call one (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat ) 

elseif (matnum.eq. 2) then 

call two (matnum, dt, nodes, epsi, k, den, cp, initT, L,T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum.eq. 3) then 

call three (matnum, dt, nodes, epsi, k,den, cp, initT, L,T, n, 
+ time, qcv, qcond, qrad, backT, heat ) 

elseif (matnum.eq. 4) then 

% call four (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat ) 

elseif (matnum.eq. 5) then 

call five (matnum, dt, nodes, epsi, k,den,cp, initT, L,T,n, 

+ time, qcv, qcond, qrad, backT, heat) 

endif 

diff (2)-T (n, 1) -templmt (1) 
root (2) -L (1) 


* While loop to call the heat equation solver. The while loop 

* reiterates until the temperature of the surface lies above the 

* frosting condition (40F), to within the specified convergence 

* criteria (error - 0.009). 


iter-2 

dowhile ( (error .gt . toler) .and. (iter . It .maxi ter) ) 
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iter - iter + 1 

root < iter) -root (iter-1) - ( (dif f (iter-1) • (root (iter-2 ) 

+ -root (iter-1 ))) / (dif f (iter-2) -dif f (iter-1 ) ) ) 

L (1) -root (iter) 

if (matnum.eq. 1) then 

call one (matnum, dt f nodes, epsi, k, den, cp, initT, L,T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum.eq. 2) then 

call two (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum.eq. 3) then 

call three (matnum, dt, nodes, epsi, k,den, cp, initT, L,T, n, 
+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum.eq. 4 ) then 

call four (matnum, dt, nodes, epsi, k, den,cp, initT, L,T,n, 

♦ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum. eq. 5) then 

call five (matnum, dt, nodes, epsi, k, den, cp, initT, L,T,n, 

+ time, qcv, qcond, qrad, backT, heat) 

endif 

dif f (iter ) -T (n, 1) -templmt (1) 

error-abs ( (root (iter) -root (iter-1) ) /root (iter) ) 
enddo 


• formatted file output statements to write temperature profiles 

• at key nodes of the TPS material stackup. 


time ( 1 ) -0 . 
do j«2,n 

time ( j ) -time { j-1 ) +dt 
enddo 

write (91,2) 
do j-l,n 

write (91,4) time { j ) , qcv ( j ) , qrad ( j } , qcond (j ) 
enddo 

if (matnum.eq. 1) then 
write (90,1) 
do j-l,n 

write (90,3) time ( j ) ,T(j,l),T(j, nodes (1 ) ) 
enddo 

elseif (matnum.eq. 2) then 
write (90,5) 
do j-l,n 

write (90,6) time <j),T(j,l),T(j, nodes ( 1 ) ) , T ( j , nodes (1) 
1+nodes (2) ) 
enddo 

elseif (matnum.eq. 3) then 
write (90, 7) 
do j«l,n 

write (90, 8) time ( j ) , T { j, 1) ,T (j , nodes (1) ), T (j, nodes (1) 
1+nodes (2) ) ,T ( j , nodes (1) +nodes (2) +nodes (3) ) 
enddo 

elseif (matnum.eq. 4) then 
write (90,9) 
do j-l,n 

write (90, 10) time (j ) , T (j , 1 ), T (j , nodes (1 ) ) ,T (j , nodes (1) 
1+nodes (2) ) , T ( j , nodes (1) +nodes (2) +nodes (3) ) , T ( j , nodes (1) +nodes (2) 
1+nodes (3) +nodes (4 ) ) 
enddo 

elseif (matnum. eq . 5) then 
write (90,11) 
do j-l,n 

write (90, 12) time (j),T(j,l),T(j, nodes ( 1 ) ) , T ( j , nodes (1) 
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l+nodes (2) ) , T ( j , nodes (1) +nodes (2) +nodes (3) ) , T ( j, nodes < 1 ) +nodes (2) 
l+nodes (3) +nodes <4) ) ,T ( j , nodes (1 ) +nodes (2) +nodes (3) +nodes (4)4 
Inodes (5) ) 
enddo 
endif 

write (12, *) Ml ) 


* Formatting parameters for file output. 


1 format (2x, 'time ( j ) * , 7x, *T ( j , 1) • , 5x, 'T ( j , nodes (1) ) * ) 

2 format (2x, 'time' ,7x, 'qconv* , 5x, 'qrad 1 , 5x, 'qcond’ ) 

3 format (f 12.2, f 12.2, f 12. 2) 

4 format (fl2 .2, 3x, f 11 .2, 3x, f 11 .2, 3x, fll .2) 

5 format (5x, ■ time* , 5x, *T ( j, 1) 1 , 5x, ' T < j, nodes 1) * , 5x 
1, *T ( j , totnodes) ') 

6 format (f 12.2, f 12.2, f 12.2, f 12. 2) 

7 format <2x, ’time(j) , ,7x,'T(j,l) \5x, 'T (j , nodes (1 ) ) \ 2x 
1, ' T ( j , suml ) ' , 5x, ' T ( j , sum2) ' ) 

8 format <f 8.2, 7x, f 8 . 2, 5x, f8 . 2, 4x, f B . 2, 5x, f8 . 2) 

9 format (2x, ' time ( j ) * , 7x, * T < j , 1 ) ',5x, *T ( j, nodes (1) ) ', 2x 
1 , ' T ( j , suml) * , 5x, 'T ( j , sum2) 1 , 5x, 'T ( j , totnodes) 1 ) 

10 format (f 8 . 2, 7x, f B . 2 , 5x, f 8 .2, 4x, f B . 2, 5x, f 8 . 2, 5x, f 8 . 2) 

11 format (2x, ’time(j) * , 7x, ' T ( j , 1) ' , 5x, *T (j , nodes (1) ) \2x 

1, 'T ( j , suml) ' , 5x, *T ( j , sum2) ’ , 5x, * T ( j , sum3) ’ , 5x, 'T (j , totnodes) ' ) 

12 format (f 8.2, 7x,fB 2 , 5x, f 8 . 2, 4x, f 8 . 2 , 5x, f 8 . 2, 5x, f 8 . 2 , 5x, f 8 . 2 ) 

end 


Pane 64 of 81 


Stephen James Izon 


Improvements in Thermal Protection Systems Sizing, CanabiUtiei 


Appendix K: ‘ground2.r FORTRAN Source Code 

program ground 


• Program to solve the one dimensional heat equation through a stackup 

• of disparate TPS materials. This program finds the optimum thickness • 

• of the surface layer of the TPS material to reduce the TPS unit • 

• weight, while eliminating frosting at the surface of the TPS stackup. • 


implicit none 

integer i, j ,n, nodes (6) , suml, sum2, sum3, totnodes,matnum, count 
integer iter f maxiter 

real* 8 dt, k (6) , den (6) , cp (6) , initT, L(6) , dx (6) , sigma, epsi,eps 
1, alpha (6) , lam (6) , a MOO) ,b (400) , c (400) , T {10000, 400) , normf 
1 , normdelT, f(400),delT(400), twonorm, time (10000) 

1 , qcond (400) , qrad (400) ,qcv (400) , g (24) , obj , x (12) , dTemp 
1, kair,Tb (1000) , error, diff (1000) ,templmt (10) ,backT,heat 
1, root (1000) , toler 


• Obtains information from the file ’inputs. in’ which contains TPS • 

• material properties, used in the solution of the heat equation. • 


open (10, file- ' inputs .in’) 
read(10,*) matnum 
read(10,*) dt 
do i-l,matnum 
readdO,*) nodes(i) 
enddo 

read(10,*) epsi 
do i-l,matnum 
readdO, •> k(i) 
enddo 

do i-l,matnum 
readdO,*) den(i) 
enddo 

do i«l,matnum 
readdO,*) cp(i) 
enddo 

readdO,*) initT 
do i-l,matnum 
readdO,*) L(i) 
enddo 

do i«l,2*matnum 
readdO,*) templmt(i) 
enddo 


1 /number of materials 
! / time step 

! /number of nodes per material 
! /emmissivity of surface 
’/thermal conductivities 

! /densities 

’/specific heats 
! /initial temperatures 
! /thicknesses 

’/temp limits for each material 


• Writes the density of the surface TPS material for use in calculating 

• the TPS material unit weight in the perl script. * 

• Also writes the heat transfer rate and the backface fuel temperature • 

• to files for use in the one-dimensional heat equation solver. 


open (40, file- 'heat' ) 
read (40, * ) heat 
close (40) 

open (41, file- ' fuel temp 1 ) 
read(41,*) backT 
close (41) 

open (9, f ile- 'mater ial_density ' ) 
write(9,*) den(2) 
close (9) 
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• Initializes optimization parameters. 


toler-0 . 0001 
error-1 . 
maxiter-100 


• Solves the heat equation for using two initial guesses for the • 

• thickness of the TPS. Values obtained for the objective function • 

• will be sent to the next while loop which implements the * 

• Newton-Raphson method to converge to the solution. * 


open ( 4 4 , f i le* 1 opt ’ } 


if (matnum.eq . 1) then 

call one (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat ) 

elseif (matnum. eq. 2) then 

call two (matnum, dt, nodes, epsi, k, den, cp, initT, L, T,n, 

+ time, qcv,qcond, qrad, backT, heat) 

elseif (matnum. eq. 3) then 

call three (matnum, dt, nodes, epsi, k, den, cp, initT, L,T, n, 
+ time, qcv, qcond, qrad, backT, heat ) 

elseif (matnum.eq. 4) then 

call four (matnum, dt, nodes, epsi, k, den, cp, initT, L,T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum.eq. 5) then 

call five (matnum, dt, nodes, epsi, k, den,cp, initT, L,T,n, 

+ time, qcv, qcond, qrad, backT, heat) 

endif 

diff (1)«T (n, l)-templmt (1) 
root (1)-L(2) 

write(44,*) diff(l) 

L (2) »0 . 006 

if (matnum.eq. 1) then 

call one (matnum, dt , nodes, epsi , k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum. eq. 2) then 

call two (matnum, dt, nodes, epsi, k, den, cp, initT,L, T,n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum. eq. 3) then 

call three (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 
+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum.eq. 4) then 

call four (matnum, dt, nodes, epsi , k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum.eq. 5) then 

call five (matnum, dt, nodes, epsi, k, den, cp, initT, L,T,n, 

♦ time, qcv, qcond, qrad, backT, heat) 

endif 

diff (2)-T(n, 1) -templmt (1) 
root (2>-L(2) 

write (44 , *) diff(2) 


• While loop to call the heat equation solver. The while loop 

• reiterates until the temperature of the surface lies above the 

• frosting condition (40F), to within the specified convergence 
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* criteria {error - 0.009). 


iter-2 

dowhile ( (error .gt . toler) .and. (iter . It .maxi ter) ) 
iter - iter + 1 

root (iter) -root (iter-1) - ( (diff (iter-1) * (root (iter-2) 

♦ -root (iter-1) ) ) / (diff (iter-2 ) -diff { iter-1 ) ) ) 

L (2)»root (iter) 
if (matnum.eq.l) then 

call one (matnum, dt, nodes, epsi , k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum. eq. 2) then 

call two (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

♦ time, qcv, qcond, qrad,backT, heat) 

elseif (matnum. eq. 3) then 

call three (matnum, dt, nodes, epsi, k,den, cp, initT, L,T,n, 
+ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum. eq. 4) then 

call four (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

♦ time, qcv, qcond, qrad, backT, heat) 

elseif (matnum. eq. 5) then 

call five (matnum, dt, nodes, epsi, k, den, cp, initT, L, T, n, 

+ time, qcv, qcond, qrad, backT, heat) 

endif 

diff (iter) -T (n, 1 ) -templmt (1 ) 

error-abs ( (root (iter) -root (iter-1) ) /root (iter) ) 
write (44, *> diff (iter) 
enddo 


close (44) 


* Formatted file output statements to write temperature profiles 

* at key nodes of the TPS material stackup. 


time (l)-0. 
do j-2,n 

time ( j ) -time ( j - 1 ) +dt 
enddo 

write (91,2) 
do j«l,n 

write <91, 4) time (j) ,qcv(j) , qrad (j ), qcond (j ) 
enddo 

if (matnum.eq.l) then 
write (90, 1) 
do j-l,n 

write ( 90, 3) time ( j ) , T ( j, 1) ,T (j , nodes (1) ) 
enddo 

elseif (matnum. eq. 2) then 
write (90, 5) 
do j-l,n 

write (90, 6) time ( j ), T ( j, 1 ),T(j, nodes (l)),T(j, nodes (1) 
1+nodes (2) ) 
enddo 

elseif (matnum. eq. 3) then 
write (90, 7 ) 
do j«l,n 

wr ite (90, 8) time ( j ), T (j,l),T(j, nodes (l)),T(j, nodes (1) 
1+nodes (2) ) ,T ( j , nodes (1) +nodes (2) ♦nodes (3) ) 
enddo 
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elseif (matnum.eq.4) then 
write (90, 9) 
do j-l,n 

write (90, 10) time (j),T(j,l),T(j, nodes (1) ) ,T ( j , nodes (1) 
1+nodes (2) ) ,T ( j , nodes (1) +nodes (2) +nodes (3) ) , T ( j , nodes (1) +nodes (2) 
1+nodes (3) +nodes (4) ) 
enddo 

elseif (matnum.eq. 5) then 
write (90, 11 ) 
do j“l,n 

write (90, 12) time ( j ) ,T ( j , 1) , T (j, nodes (1) ) ,T ( j , nodes (1) 
1+nodes (2 ) ) , T ( j , nodes (1) +nodes (2) +nodes (3) ) , T ( j , nodes ( 1 ) +nodes (2) 
1+nodes (3) +nodes (4) ) ,T (j, nodes (1 ) +nodes (2) +nodes (3) +nodes (4) + 
Inodes (5) ) 
enddo 
endif 

write (12, * ) L (2) 


* Formatting parameters for file output. 


1 format <2x, ’time(j) \7x, 'T(j,D ’ ,5x, *T(j, nodes (U) ' ) 

2 format (2x, ’time* , 7x, ’qconv' , 5x, 'qrad' , 5x, ' qcond* ) 

3 format (f 12.2, f 12.2, f 12. 2) 

4 format (f 12.2, 3x,f 11. 2, 3x, f 11 .2, 3x, f 11 .2) 

5 format (5x, * time' , 5x, 'T ( j, 1 ) ',5x, ’T(j, nodesl) ’ , 5x 
1, 'T ( j , totnodes) * ) 

6 format (f 12 . 2 , f 12 . 2, f 12 . 2, f 12 .2) 

7 format (2x, 'time ( j ) * , 7x, *T ( j , 1) ' , 5x, 'T ( j , nodes (1) ) * , 2x 
1, ’ T ( j , suml ) * , 5x, *T(j,sum2) * ) 

8 format (f 8.2, 7x, f 6 . 2 , 5x, f8 . 2, 4x, f 8 . 2 , 5x, f 8 . 2) 

9 format (2x, ' time (j) ’,7x, ' T < j, 1) *,5x, *T(j, nodes (1) ) ’ , 2x 
1, 'T ( j , suml) ' , 5x, ' T ( j , sum2 ) ' , 5x, 'T ( j , totnodes) ' ) 

10 format (f 8 . 2, 7x, f 8 . 2 , 5x, f 8 . 2, 4x, f 8 . 2 , 5x, f 8 . 2 , 5x, f 8 . 2) 

11 format (2x, 'time (j ) ' , 7x, 'T ( j , 1) ' , 5x, 'T ( j , nodes (1 ) ) ’ , 2x 

1, 'T ( j, suml ) ’ , 5x, * T ( j , sum2) ' , 5x, ' T ( j , sum3) * , 5x, ' T ( j , totnodes) 1 ) 

12 format (f8.2,7x,f8.2,5x,f8.2,4x,f8.2,5x,f8.2,5x,f8.2,5x,f8.2) 

end 
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Appendix L: ‘one.r FORTRAN Subroutine 


subroutine one (matnum, dt, nodes, epsi, k, den, cp, initT, L,T,n, 
+ time, qcv, qcond, qrad, backT, heat) 


• Subroutine to solve the one dimensional heat equation for a TPS • 

• material stackup consisting of one materials. 

• This subroutine incorporates an isothermal boundary condtion at the 

• backface of the TPS material stackup (generally the cryogenic fuel 

• temperature) and also a constant heat flux boundary condition at the 

• surface of the TPS material stackup. The constant heat flux is • 

• determined solely by the temperature of the ambient air. • 


implicit none 

integer i, j , n, nodes (6) , suml , sum2, sum3, totnodes, matnum, count 
real *8 dt, k (6) , den ( 6) , cp (6) , initT, L (6) , dx (6) , sigma, epsi, eps 
1, alpha (6) ,1am (6) ,a(4 00) , b (400) , c (400) , T (10000, 400) , normf 
1, normdelT, f(400) ,delT(400), twonorm, Tradeq (10000) , time (10000) 
1, qcond (10000) , qrad (10000) , qcv (10000) , g (24 ) , obj , x < 12 ) , Tb ( 400) 
1, backT, heat 


• Calls linint subroutine to linearly interpolate convective heat flux 

• and Radiation Equilibrium temperature values using user specified • 

• time and time step for analysis. * 


call linint (n, dt, time, qcv, Tradeq) 


* Define the mathematical constants to be used in the calculations. • 


eps-l . e-6 
sigma-5 . 67061e-8 
do i-1, matnum 

dx (i) -L (i) / (nodes (i) -1 . ) 
alpha (i)-k(i)/ (den (i ) *cp(i) *1000.) 
lam (i) -alpha (i) *dt/ (dx (i) **2) 
enddo 

suml -nodes (1 ) +nodes (2 ) 
sum2-suml+nodes (3) 
sum3-sum2+nodes (4) 

totnodes -nodes (1 ) +nodes (2 ) + nodes (3) +nodes (4 ) +nodes (5) 


* Specify the initial conditions of the problem. 


do i-1, totnodes 
T (1 , i ) -initT 
enddo 

do i-1, totnodes 
Tb (i) -T (1, i) 
enddo 

Tb (totnodes) -backT 

do i-1, totnodes 
T ( 1 , i ) -Tb (i ) 
enddo 

********************«******** 4 ************+********+*******»************* 

* Initialize time array. * 

************************************************************************* 


time (1 )-0. 


Stephen James Izon 


Pace 69 of SI 


Improvements in Thermal Protection Systems Sizirte Capabilities for TCA T 


• Do loop to solve for the temperature profiles for each time step j 
************************************************************************** 


do j«2,n 


• Setup the tridiagonal matrix to be used in the Thomas Algorithm. 


a (1)— 2.*lam(l) 
b ( 1) *1 . +2 . *lam (1) 
c (1) *0 . 

do i-2, nodes (1) -1 
a<i) — l.*lam(l) 
b(i) -1 .+2. *lam (1) 
c (i) •- 1 . *lam(l) 
enddo 

a (totnodes) *0 . 
b (totnodes ) -1 . 
c (totnodes) -0 . 

time ( j ) -time < j-1) +dt 

Tb ( 1 ) «Tb ( 1 ) +2 . * lam ( 1 ) *heat *dx < 1 ) /k ( 1 ) 

call tridag (a, b, c, Tb r de IT , totnodes) 

do i-1 , totnodes 
T(j, i)-delT(i) 

Tb { i) -delT (i) 
enddo 

enddo 


• Calculation of the radiative and conductive heat f luxes, - 


do j*l,n 

qrad ( j ) «epsi*sigma*T ( j , 1 ) **A 
qcond{ j )«qcv( j) -qrad(j) 
enddo 


• »***************************’********•****************•**'******••******** 

• End of subroutine one. • 


return 

end 
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Appendix M: ‘two.f FORTRAN Subroutine 

subroutine two (matnum, dt, nodes, epsi, k, den, cp, initT, L,T,n, 
+ time, qcv, qcond, qrad, backT, heat ) 


• Subroutine to solve the one dimensional heat equation for a TPS * 

• material stackup consisting of two dissparate materials. * 

• This subroutine incorporates an isothermal boundary condtion at the 

• backface of the TPS material stackup (generally the cryogenic fuel 

• temperature) and also a constant heat flux boundary condition at the 

• surface of the TPS material stackup. The constant heat flux is • 

• determined solely by the temperature of the ambient air. • 


implicit none 

integer i, j,n, nodes (6) , suml, sum2,sum3, tot nodes, matnum, count 
real *8 dt,k(6) ,den(6) , cp(6) , initT, L(6) ,dx(6) , sigma, epsi, eps 
1, alpha (6) , lam (6) , a <4 00) ,b (400) , c (400) ,T (10000, 400) , normf 
1, normdelT, f (400) , delT (400) , twonorm, Tradeq (10000) , time (10000) 
l,qcond (10000) , qrad (10000) ,qcv (10000) , g (24 ) , obj , x (12) , Tb (400) 
l,backT,heat 


• Calls linint subroutine to linearly interpolate convective heat flux 

• and Radiation Equilibrium temperature values using user specified * 

• time and time step for analysis. * 


call linint (n, dt, time, qcv, Tradeq) 


• Define the mathematical constants to be used in the calculations. * 


eps-1 .e-6 
sigma-5 . 67061e«^8 
do i-1, matnum 

dx (i) -L (i) / (nodes (i) -1 . ) 
alpha (i)-k(i) / (den (i) *cp (i) *1000.) 
lam (i) -alpha (i) *dt/ (dx(i) **2) 
enddo 

suml-nodes (1 ) +nodes (2) 
sum2-suml+nodes (3) 
sum3-sum2+nodes (4) 

totnodes-nodes (1) +nodes (2) +nodes (3) +nodes (4 ) +nodes (5) 


• Specify the initial conditions of the problem. 


do i-l,totnodes 
T (1, i) -initT 
enddo 

do i«l,totnodes 
Tb ( i ) -T ( 1 , i ) 
enddo 

Tb (totnodes) -backT 

do i-1, totnodes 
T(l,i)«Tb<i) 
enddo 


• Initialize time array. 


time (1 ) -0 . 
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♦♦a***********-*********************************************************** 
* Do loop to solve for the temperature profiles for each time step j 


do j-2,n 


* Setup the tridiagonal matrix to be used in the Thomas Algorithm. 


a (l)--2,*lam(l) 

b(l)-l.+2.*lam(l) 

c<l)-0. 

do i-2, nodes (1) -1 
a (i) — 1 .*lam (1> 
b (i) -1 .+2 . *lam(l) 
c (i)— l.*lam(l) 
enddo 

a (nodes (1) ) —2 . *lam (2) 

b (nodes (1) > -1 .+2 . *lam (2) +2 . *lam (2) *k(l>/k(2) 
c (nodes (1) ) — 2 . *lam (2) *k (1) /k (2) 
do i-nodes (1) +1 , suml-1 
a (i) — 1 . *lam (2) 
b(i)«l.+2.*lam(2> 
c (i ) — 1 . *lam (2 ) 
enddo 

a (totnodes) *0 . 
b (totnodes) -1 . 
c (totnodes ) -0 . 

time ( j ) -time ( j-1) +dt 

Tb(l)-Tb(l) +2.*lam(l) *heat*dx(l) /k(l) 

call tridag (a, b, c,Tb, delT, totnodes) 

do i-1, totnodes 
T ( j , i) -delT (i) 

Tb (i ) -delT (i) 
enddo 

enddo 


* Calculation of the radiative and conductive heat fluxes. 


do j-l,n 

qrad { j ) -epsi*sigma*T (j, 1) **4 
qcond { j ) -qcv (j)-qrad(j) 
enddo 


* End of subroutine two. 


return 

end 
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Appendix N: ‘three.f FORTRAN Subroutine 

subroutine three (matnum, dt, nodes, epsi, k, den, cp, initT, L,T, n, 
+ time, qcv, qcond, qrad, backT, heat) 


• Subroutine to solve the one dimensional heat equation for a TPS • 

• material stackup consisting of three dissparate materials. * 

• This subroutine incorporates an isothermal boundary condtion at the * 

• backface of the TPS material stackup (generally the cryogenic fuel • 

• temperature) and also a constant heat flux boundary condition at the • 

• surface of the TPS material stackup. The constant heat flux is • 

• determined solely by the temperature of the ambient air. • 


implicit none 

integer i, j , n, nodes (€) , suml, sum2, sum3, totnodes, matnum, count 
real *8 dt, k (6) , den (6) , cp(6) , initT, L (6) , dx (6) , sigma, epsi, eps 
1, alpha (6) , lam (6) , a (400) ,b (400) , c (400) ,T (10000,400) ,normf 
1, normdelT, f (400) , delT (400) , twonorm, Tradeq (10000) , time (10000) 
1, qcond (10000 ), qrad ( 10000) , qcv (10000) ,g (24) ,obj,x (12) ,Tb (400) 
1, backT, heat 


• Calls linint subroutine to linearly interpolate convective heat flux 

• and Radiation Equilibrium temperature values using user specified • 

• time and time step for analysis. * 


call linint (n, dt, time, qcv, Tradeq) 


* Define the mathematical constants to be used in the calculations, * 


eps-1 .e-6 
sigma-5 . 67061e-8 
do i-1, matnum 

dx (i) -L (i) / (nodes (i) -1 . ) 
alpha (i) -k (i) / (den (i) *cp (i) *1000. ) 
lam (i) -alpha (i) *dt/ (dx (i) **2) 
enddo 

suml -nodes (1 ) +nodes (2) 
s um2 -suml ♦ nodes (3) 
sum3-sum2+nodes (4) 

totnodes-nodes (1 ) +nodes (2) +nodes (3) +nodes (4 ) +nodes (5) 


* Specify the initial conditions of the problem. 


do i-1, totnodes 
T (1 , i) -initT 
enddo 

do i-1, totnodes 
Tb (i) -T (1, i) 
enddo 

Tb (totnodes) -backT 

do i-1, totnodes 
T (1, i) -Tb (i) 
enddo 


* Initialize time array. 


time (1 ) -0 . 
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* Do loop to solve for the temperature profiles for each time step j 


do j-2,n 


* Setup the tridiagonal matrix to be used in the Thomas Algorithm. 


a(l> — 2.*lam(l) 
b(l)-l.+2.*lam<l) 
c (1) -0 . 

do i-2, nodes < X ) -1 
a (i) — 1 .*lam (1) 
b(i)-l.+2.*lam(l) 
c (i) — 1 . *lam (1) 
enddo 

a (nodes (1) ) --2 . *lam (2) 

b (nodes (1) > -1 . +2 . *lam (2 ) +2 . *lam (2) *k(l)/k(2> 
c (nodes (1) ) — 2 . *lam (2) *k(l> /k(2) 
do i-nodes (1) +1, suml-1 
a (i) —1 . *lam (2) 
b (i) -1 . +2 . *lam (2) 
c (i) — 1 . *lam (2) 
enddo 

a (suml) — 2 . *lam (3) 

b (suml ) -1 . +2 . *lam (3) +2 . Mam (3) *k(2)/k(3) 
c (suml)— 2. *1 am (3) *k(2) /k(3) 
do i-suml+1 , sum2-l 
a fi) —1 . *lam (3) 
b(i)«l.+2.*lam(3) 
c (i)— 1 . *lam (3) 
enddo 

a (totnodes) -0 . 
b (totnodes) -1 . 
c (totnodes) -0 . 

time ( j ) -time (j-1) +dt 

Tb (l)-Tb(l) +2. *lam(l) *heat*dx(l) /k(l) 

call tridag (a, b, c,Tb, del T, totnodes) 

do i-1 , totnodes 
T(j, i)-delT(i) 

Tb (i) -delT (i ) 
enddo 

enddo 


* Calculation of the radiative and conductive heat fluxes. 


do j-l,n 

qrad ( j ) -epsi* sigma*T < j , 1) **4 
qcond ( j ) -qcv ( j ) -qrad ( j ) 
enddo 


* End of subroutine three. 


return 

end 
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Appendix O: ‘four.r FORTRAN Subroutine 

subroutine four (matnum, dt, nodes , eps i , k, den, cp, initT,L, T, n, 
+ time, qcv, qcond, qrad, backT, heat) 


• Subroutine to solve the one dimensional heat equation for a TPS • 

• material stackup consisting of four dissparate materials. * 

• This subroutine incorporates an isothermal boundary condtion at the 

• backface of the TPS material stackup (generally the cryogenic fuel 

• temperature) and also a constant heat flux boundary condition at the 

• surface of the TPS material stackup. The constant heat flux is * 

• determined solely by the temperature of the ambient air. * 


implicit none 

integer i, j, n, nodes (6) , suml , sum2 f sum3, totnodes,matnum, count 
real *8 dt, k(€) ,den(6) ,cp(6) ,initT,L(6) ,dx(6) , sigma, eps i, eps 
1, alpha (6) ,1am (6) , a (400) ,b (400) , c (400) , T (10000, 400) r normf 
1, normdelT, f (400) , delT (400) , twonorm, Tradeq (1 0000) , time (10000) 
l,qcond (10000) , qrad (10000) , qcv (10000) ,g (24) , ob j , x (12) , Tb (400) 
1, backT, heat 


• Calls linint subroutine to linearly interpolate convective heat flux 

• and Radiation Equilibrium temperature values using user specified • 

• time and time step for analysis. • 


call linint (n, dt, time, qcv, Tradeq) 


* Define the mathematical constants to be used in the calculations. * 


eps-l.e-6 
sigma-5 . 67061e-8 
do i*l, matnum 

dx (i) *L (i) / (nodes (i) -1 . ) 
alpha (i) -k (i) /(den (i) *cp(i)* 1000.) 
lam (i) -alpha (i) *dt/ (dx (i) **2) 
enddo 

suml-nodes (1) Anodes (2) 
sum2-suml+nodes (3) 
sum3*sum2+nodes (4) 

totnodes-nodes (l)+nodes (2)+nodes (3) +nodes (4) +nodes (5) 


• Specify the initial conditions of the problem. 


do i«l,totnodes 
T (1 , i ) -initT 
enddo 

do i«l,totnodes 
Tb(i)-T(l,i) 
enddo 

Tb (tot nodes) -backT 

do i-l,totnodes 
T(l,i)-Tb(i) 
enddo 


* Initialize time array. 


Stephen James Izon 


Pane 75 of 81 


Improvements in Thermal Protection Systems Sizins Capabilities for TCAT 


time (1) -0 . 


* Do loop to solve for the temperature profiles for each time step j 


do j-2,n 


• Setup the tridiagonal matrix to be used in the Thomas Algorithm. 


a (D— 2.*lam(l) 
b(l)-l.+2.*lam(l) 
c (1) -0. 

do i-2, nodes (1) -1 
a (i)--l . *lam(l) 
b(i)-l .+2. *lam(l) 
c (i) — 1 . *lam (1) 
enddo 

a (nodes (1 ) ) — 2 . *lam (2) 

b (nodes (1 ) ) -1 . +2 . *lam (2) +2 . *lam (2) *k(l) /k(2) 
c (nodes (1 ) ) — 2 . *lam (2) •k(l)/k(2) 
do i-nodes (1) +1 , suml-1 
a ( i ) — — 1 . *lam (2) 
b (i ) -1 . +2 . *lam (2) 
c (i) --3 . *lam (2) 
enddo 

a (suml) — ? . *lam (3) 

fc (sum! ) . +? . *lam (3> +2 . *lam (3) *k(2)/k(3) 
c (suml)— 2. *lam(3) *k(2)/k(3) 
do i-suml+l,sum2-l 
a(i> — l.*lam{3) 
b (i ) -1 . +2 . *lam (3) 
c(i) — l.*lam(3) 
enddo 

a (sum2) — 2 . *lam (4) 

b(sum2>«l.+2,*lam(4)+2.*lam<4)*k(3) /k{4) 
c (sum2)— 2.*lam(4) *k<3)/k(4> 
do i-sum2+l , sum3-l 
a(i)— l.*lam(4) 
b <i) -1 . +2 . *lam(4) 
c(i)— l.*lam(4) 
enddo 
a (sum3) - 0 . 
b (sum3) -1 . 
c (sum3) -0 . 

time { j ) -time ( j-1) +dt 

Tb (1) -Tb (1) +2 . *lam (1 ) *heat*dx (1 ) /k (1) 

call tridag (a, b, c, Tb, delT, totnodes) 

do i-1, totnodes 
T ( j , i) -delT (i) 

Tb(i)-delTU) 

enddo 

enddo 


* Calculation of the radiative and conductive heat fluxes. 


do j-1, n 

qrad ( j ) -epsi*sigma*T (j , 1) **4 
qcond ( j ) -qcv( j ) -qrad { j ) 
enddo 
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* End of subroutine four. 


return 

end 
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Appendix P: ‘flve.f FORTRAN Subroutine 

subroutine five (matnum, dt , nodes, epsi, k, den, cp, initT, L, T,n, 
+ time, qcv, qcond, qrad, backT, heat) 


• Subroutine to solve the one dimensional heat equation for a TPS • 

• material stackup consisting of five dissparate materials. * 

• This subroutine incorporates an isothermal boundary condtion at the 

• backface of the TPS material stackup (generally the cryogenic fuel 

• temperature) and also a constant heat flux boundary condition at the 

• surface of the TPS material stackup. The constant heat flux is • 

• determined solely by the temperature of the ambient air. • 


implicit none 

integer i, j,n, nodes (6) , suml,sum2,sum3, totnodes, matnum, count 
real *8 dt, k(6) ,den (6) , cp (€) , initT, L (6) ,dx (6) , sigma, epsi, eps 
1, alpha (6) , lam (6) , a (400) ,b (400) , c (400) ,T (10000, 400) ,normf 
l,normdelT, f (400) ,delT(400) , twonorm, Tradeq (10000) , time(lOOOO) 
1, qcond (10000) , qrad (10000) , qcv (10000) ,g (24) , ob j,x (12) ,Tb (400) 
1, backT, heat 


• Calls linint subroutine to linearly interpolate convective heat flux 

• and Radiation Equilibrium temperature values using user specified • 

• time and time step for analysis. • 


call linint (n, dt, time, qcv, Tradeq) 


• Define the mathematical constants to be used in the calculations. • 


eps-l . e-6 
sigma-5. 67061e-8 
do i-1, matnum 

dx (i) -L (i) / (nodes (i) -1 . ) 
alpha (i) *k (i) / (den (i)*cp(i)*1000.) 
lam (i) -alpha (i) *dt/ (dx (i) ** 2 ) 
enddo 

suml-nodes (1 ) +nodes (2 ) 
sum2-suml+nodes (3) 
sum3-sum2+nodes (4) 

totnodes-nodes (1) +nodes (2) +nodes (3)+nodes (4)+nodes (5) 


• Specify the initial conditions of the problem. 


do i-1, totnodes 
T (l r i) -initT 
enddo 

do i-1, totnodes 
Tb (i ) -T (1 , i ) 
enddo 

Tb (totnodes) -backT 

do i-1, totnodes 
T ( 1, i ) -Tb (i ) 
enddo 


• Initialize time array. 


time (l)-0. 
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# Do loop to solve for the temperature profiles for each time step j 


do j-2,n 


* Setup the tridiagonal matrix to be used in the Thomas Algorithm. 


a (1)— 2 . *lam (1) 
b(l)«l.+2.*lam(l) 
c (1) -0 . 

do i«2, nodes (1) -1 
a (i) — 1 . *lam (1) 
b(i)-l.+2.*lam(l) 
c (i) — 1 . *lam(l) 
enddo 

a (nodes (1) ) — 2. *1 am(2) 

b (nodes (1) ) -1 .+2 .*lam (2) +2 . *lam (2) *k(l)/k<2) 
c (nodes (1) ) — 2 . *lam (2) *k (1 ) /k (2) 
do i-nodes <1)+1, suml-1 
a (i) — 1 . *lam (2) 
b (i ) -1 . +2 . *lam (2) 
c(i) — l.*lam<2) 
enddo 

a (suml) —2 . *lam{3) 

b (suml) -1 , +2 . *lam(3) +2 . *lam (3) *k (2) /k(3) 
c (suml) —2 . *lam (3) *k (2) / k (3) 
do i-suml+1, sum2-l 
a (i) — 1 . *lam (3) 
b (i) -1 .+2 . *lam(3) 
c(i>— l.*lam(3) 
enddo 

a (sum2) — 2 . *lam (4 ) 

b(sum2)-l.+2.*lam(4)+2.*lam<4) *k(3)/k(4) 
c (sum2) —2 . *lam (4) *k (3) /k (4) 
do i-sum2+l r sum3-l 
a (i) *-l . *lam (4) 
b(i)*l.+2. *lam(4) 
c (i) —1 . *lam (4) 
enddo 

a { sum3 ) --2 . +lam (5) 

b (sum3) «1 . +2 . *lam{5) +2 . *lam (5) *k (4) /k{5> 
c (sum3)*-2. *lam(5)*k(4)/k<5) 
do i-sum3+l , totnodes-1 
a ( i) . *lam (5) 
b(i)-l .+2. *lam (5) 
c (i) — 1 . *lam (5) 
enddo 

a (totnodes) -0 . 
b (totnodes) *1 . 
c (totnodes) -0 . 

time { j ) -time ( j -1) +dt 

Tb ( 1 ) «Tb ( 1 ) +2 . *lam ( 1 ) *heat*dx (1) /k (1) 

call tridag (a,b, c, Tb , delT, totnodes) 

do i-1, totnodes 
T{ j, i) -delT (i) 

Tb (i)-delT(i) 
enddo 

enddo 


% Calculation of the radiative and conductive heat fluxes. 
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do j«l,n 

qrad ( j ) -epsi*sigma*T { j , 1) **4 
qcond ( j ) *qcv { j ) -qrad (j) 
enddo 


* End of subroutine four. 


return 

end 
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